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ABSTRACT 

Accurately understanding the interior structure of extra-solar planets is critical for inferring their 
formation and evolution. The internal density distribution of a planet has a direct effect on the star- 
planet orbit through the gravitational quadrupolc field created by the rotational and tidal bulges. 
These quadrupoles induce apsidal precession that is proportional to the planetary Love number (fop, 
twice the apsidal motion constant), a bulk physical characteristic of the planet that depends on the 
internal density distribution, including the presence or absence of a massive solid core. We find that 
the quadrupole of the planetary tidal bulge is the dominant source of apsidal precession for very hot 
Jupiters (a < 0.025 AU), exceeding the effects of general relativity and the stellar quadrupole by more 
than an order of magnitude. For the shortest-period planets, the planetary interior induces precession 
of a few degrees per year. By investigating the full photometric signal of apsidal precession, we find 
that changes in transit shapes are much more important than transit timing variations. With its long 
baseline of ultra-precise photometry, the space-based Kepler mission can realistically detect apsidal 
precession with the accuracy necessary to infer the presence or absence of a massive core in very hot 
Jupiters with orbital eccentricities as low as e ~ 0.003. The signal due to fo p creates unique transit 
light curve variations that are generally not degenerate with other parameters or phenomena. We 
discuss the plausibility of measuring fo p in an effort to directly constrain the interior properties of 
extra-solar planets. 
Subject headings: stars: planetary systems 
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1. INTRODUCTION 

Whether studying planets within our solar system or 
planets orbiting other stars, understanding planetary in- 
teriors represents our best strategy for determining their 
bulk composition, internal dynamics, and formation his- 
tories. For our closest neighbors, we have had the luxury 
of sending spacecraft to accurately measure the higher- 
order gravity fields of these objects, yielding invaluable 
constraints on their interior density distributions. Using 
these observations, we have been able, for instance, to in- 
fer the presence of large cores, providing s upport for th e 
core-accretion theory of planet formation (|Guillotll2005l) . 
Study of planets outside our solar system, however, has 
necessitated the development and usage of more indi- 
rect techniques. Nevertheless, as the number of well- 
characterized extra-solar planets grows, we gain more 
clues that help us answer the most fundamental ques- 
tions about how planets form and evolve. 

Guided by our current understanding of planetary 
physics, we have begun to study the interiors of extra- 
solar planets. This endeavor has been dominated by a 
model-based approach, in which the mass and radius of 
a planet are measured using radial velocity and transit 
photometry observations, and the interior properties are 
inferred by finding the model most consistent with those 
two observations. This strategy clearly requires a set of 
assumptions, not the least of which is that the physical 
processes at work in extra-solar planets are just like those 
that we understand for our own giant planets. While it 
does seem that this approach is adequate for explaining 
most of the known transiting planets, there does exist 
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a group of planets for which the usual set of assump- 
tions are not capable of reproducing the o bservations 
(e.g.. lGuillot et alj|2006t IBurrows et alj|2007l ). These are 
the planets with so-called positive "radius anomalies" , in- 
cluding the first-discovered transiting planet HD 209458b 
(jCharbonneau et al.l 120001 ). Though most of these plan- 
ets can be explained by adjusting different pieces of the 
interior physics in the models (including opacities, equa- 
tions of state, and heat deposition), it is currently impos- 
sible to discern which combination of these possible ex- 
planations is actual ly responsible for their observed sizes 
ijGuillot et afll2006l) . 

Additional uncertainties also exist for planets at the 
other end of the size spectrum. For the group of under- 
sized extra-solar planets, such as HD 149026b, the canon- 
ical approach is to give the planet a massive highly con- 
densed core of heavy elements in order to match the ob- 
served radius. This approach also provides a first order 
estimate of the planet's bulk composition, in terms of 
its fraction of heavy elements. There is also the added 
complication of how the assumed state of differentiation 
affects the inferred composition and predicted structure 
(jBaraffe et al.ll2008h . 

Currently, the most promising approach to modeling 
the distinctive features of extra-solar planet interiors is 
to study the known transiting planets as an ensemble. 
The group can be used to develop either a single con- 
sistent model that reproduces all the observations (e.g., 
iGuillot et al.| [2~006) or to showcase the po ssible diversity 
in mo del parameters (e.g., opacities, as in IBurrows et al.l 
2007). Surely, a model-independent measure of interior 
structure would be valuable in order to begin disentan- 
gling otherwise unconstrained physics. 

The idea of obtaining direct structural measurements 
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for distant objects is by no means a new one. For 
decades, the interiors of eclipsing binary stars have been 
measured by observing "apsidal motion," i.e. precession 
of the orbit due to the non-point-mass component of the 
gravit ational field (|Russel]||192l ICowlindll93l ISternd 
H939allbT ). The signal of the changing orbit is encoded in 
the light curves of these systems by altering the timing of 
the primary and secondary eclipses. From these eclipse 
times, it is straightforward to determine the so-called ap- 
sidal motion constant which then constrains the allowed 
interior density distributions. Interior measurements in- 
ferred from apsidal precession were among the first indi- 
cations that stars were highly centrally condensed. While 
it seems non-intuitive, we show in this paper that we can 
use a similar technique to measure the interior proper- 
ties of very hot Jupiters. Most surprisingly, the interior 
structure signal for very hot Jupiters actually dominates 
over the signal from the star, yielding an unambiguous 
determination of planetary interior properties. 

Our theoretical analysis is also extended to full sim- 
ulated photometry in order to explore the observability 
of apsidal precession. We show that this precession is 
observable by measuring the subtle variations in tran- 
sit light curves. The photometric analysis is focused on 
the data expected from NASA's Kepler mission, which 
succe s sfully launched o n March 6, 2009 (|Borucki et al.1 
120031 iKochet all [20061. Kepler will obtain exquisite 
photometry on ^100,000 stars, of which about 30 are 
expecte d to host hot Jupiter s with periods less than 
3 days (|Beattv fc Gaudill2008t ). Kepler has the poten- 
tial to measure the gravitational quadrupoles of very hot 
Jupiters though the technique described below. If suc- 
cessful, this will constitute a major step towards an un- 
derstanding of the diversity of planetary interiors. 

In Section 2, we describe the background theory that 
connects interior structure and orbital dynamics and ex- 
plore which effects are most important. Section 3 ap- 
plies this theory to the observable changes in the transit 
photometry, including full Kepler simulated light curves. 
We show in Section 4 that the signal due to the planetary 
interior has a unique signature. Other methods for infer- 
ring planetary interior properties are discussed in Section 
5. The final section discusses the important conclusions 
of our work. 

2. BACKGROUND THEORY 

2.1. Coordinate System and Notation 

The internal structure of very hot Jupiters can be 
determined by observing changes in the planet's orbit. 
These changes can be described in terms of two general 
types of precession. Apsidal precession refers to rota- 
tion of the orbital ellipse within the plane of the orbit. 
It is characterized by circulation of the line of apsides, 
which lies along the major axis of the orbit. Nodal pre- 
cession, on the other hand, occurs out of the plane of 
the orbit and refers to the orbit normal precessing about 
the total angular momentum vector of the system. For 
typical very hot Jupiter systems with no other planets, 
apsidal precession has a much stronger observable signal 
than nodal precession (see Section HTTj) . so we focus our 
discussion on the simpler case of a fixed orbital plane. 

As is typical for non-Keplerian orbits, the star-planet 
orbit is described using osculating orbital elements that 
change in time. We identify the plane of the sky as the 



reference plane and orient the coordinate axes in the 
usual way such that the sky lies in the x-z plane with 
the y-axis pointing at Earth. The intersection of the or- 
bital plane and the reference plane is called the line of 
nodes, but without directly resolving the system, there is 
no way to determine the orientation of the line of nodes 
with respect to astronomical North; thus, the longitude 
of the ascending node, Q, cannot be determined. Given 
this degeneracy, we simplify the description by orienting 
the z-axis to lie within the plane spanned by the orbit 
normal and the line-of-sight. The angle between the line 
of sight and the orbit normal is i, the inclination. The 
x-axis is in the plane of the sky and is the reference line 
from which the argument of periapse (lo) is measured (in 
the standard counter-clockwise sense) . For this choice of 
coordinates, the argument of periapse and longitude of 
periapse (w) are equivalent. Given this coordinate sys- 
tem, transit centers occur when the planet crosses the y-z 
plane; this point lies 90° past the reference x-axis, and 
thus primary transits occur when the true anomaly, /, 
satisfies ftr + ^tr = 90°, where the subscript tr indicates 
the value at transit center. 1 

Throughout this paper, we refer to parameters of the 
star (mass, radius, etc.) with subscripts of "*" and 
parameters of the planet with subscripts of "p" . For 
evaluation of various equations, we will take as fidu- 
cial values the mass ratio M p /M* = 10 -3 , the radius 
ratio R p /R„, =0.1 (though some low density planets 
have radius ratios greater than 1/6), and the semi-major 
axis in stellar radii a/R* = 6, typical for very hot 
Jupiters, which we define as planets with semi-major axes 
a < 0.025 AU (see Table 1) 2 In this definition, we de- 
viate from iBeattv fe Gaudil (|2008l ). who define very hot 
Jupiters as planets with periods less than 3 days. These 
authors estimate that Kepler will find ~30 such planets, 
of which ~16 will be brighter than V=14 (T. Beatty, 
pers. comm.). Since our definition is more stringent, our 
technique will be applicable to fewer Kepler planets. 

2.2. Rotational and Tidal Potentials 

It is well known from classical mechanics, that if stars 
and planets are considered to be purely spherical masses, 
then they will obey a simple r~ 2 force law and hence ex- 
ecute closed elliptical orbits. Non-spherical mass effects 
are caused by the application of external potential(s): 
the centrifugal potential of spinning bodies causes rota- 
tional flattening and the tidal potential of a nearby mass 
raises tidal bulges. Rotational and tidal bulges create 
gravitational quadrupolc fields (r~ 3 ) that lead to orbital 
precession. 

The complex subject of how planets 3 respond to ap- 
plied potenti als is encapsulated in th e so-called theory 
of figures (Z harkov fe Trubitsvnlll97c? h As long as the 
distortions are small, we can simplify the problem by ig- 
noring the small interaction terms between the tidal and 
rotational potentials; in this paper, we thus restrict our- 
selves to the first order theory, where the two planetary 

1 In elliptical orbits, if the inclination is not 90°, the photometric 
mi nima do not e xactly coincide wit h the planetary conjunctions. 
Sec Kopal (1959), p. 388 and section T3. 31 below. 

2 Throughout this work, we do not distinguish between Mtot 
and M* , since M p < M» . 

3 For clarity, in these sections we focus on the planetary shape, 
though the derivations are also valid for stars. 
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responses simply add. Even in the linear case, the way 
the fluid planet responds depends on the full radial den- 
sity structure of the planet. The planetary response is 
conveniently captured in a single variable k 2pi using the 
definition 

\q ad (R p ) = k 2p V»™(R p ) (1) 

where k 2p is the Love number of the planet, which is 
just a constant of proportionality between the applied 
second degree potential field V^ w and the resulting field 
that it induces V r 2 lnd at the surface of the planet. Due to 
the orthogonality of the Legendre polynomials used to 
express the gravity field, if the planet is responding to a 
second degree harmonic field, then only the second degree 
harmonic of the planet's gravity field is altered, to first- 
order. Thus, k 2p is a measure of how the redistribution 
of mass caused by the applied potential actually affects 
the external gravity field of the planet. In the stellar 
literature, the symbol k 2 is used for the apsidal motion 
constant, which is half of the secul ar /fluid Love n umber 
that we use throughout this paper (|Sterndll939af ). 

The Love number k 2 is an extremely useful parameter- 
ization, as it hides the complex interactions of a planet 
and an applied potential in just a single number. The 
process of calculating k 2 of a fluid object (like stars 
and gas giants), from the interior density distribution is 
fairly straightforward and outlined in several places (e.g., 
Sterne 1939a; Kopal 1959). Objects with most of their 
mass near their cores, like stars, have very low k 2 values 
(^0.03 for main sequence solar-like stars, IClaretl [l995h 
since the distorted outer envelope has little mass and 
therefore little effect on the gravity field. Planets have 
much flatter density distributions, and thus distortions 
of their relatively more massive outer envelopes greatly 
affect the gravity field. At the upper extreme lies a uni- 
form density sphere, which has k 2 = 3/2. In this way, 
k 2 can be thought of as a measure of the level of central 
condensation of an object, with stronger central conden- 
sation corresponding to smaller k 2 . 

By examining the variations in k 2 for giant planets 
within our own Solar System, we can gain a feel for its 
expected values and how sensitive it is to internal struc- 
ture. The n = l polytrope is commonly used to approx- 
imate the density structure of (cold) gas giant planets; 
it has k 2 k, 0.52 (|Kopall fl959f h This can be compared 
to the value determined from the gravity measurements 
of Jupiter, where k 2 j ~ 0.49. Even though Jupiter may 
have a 10 Earth mass core, it is small in comparison to 
Jupiter's total mass, and thus it has minor effect on the 
value of k 2 . Saturn, on the other hand has a roughly 20 
Earth mass core and is less than 1/3 of Jupiter's mass. As 
a result, the presence of Saturn's core is easily seen in the 
value of its Love number k 2 s ~ 0.32. From this, we can 
see that planets with and without significant cores differ 
in k 2v by about ~ 0.1. This can also be inferred from 
iBarnes fc Fortnevl (|2003h by using the Darwin-Radau re- 
lation to c onvert the moment of inert ia factor to k 2 . Fur- 
thermore, iBodenheimer et all (|2001h list the moment of 
inertia factors of various planet models of HD 209458 b 
and r Bootis b, which correspond to a range of k 2p values 
from ~0.1 to ~0.6. 

Current methods for inferring the internal structures 
of extra-solar planets combine measurements of the mass 
and radius with a model to obtain estimates of the 



planet's implied composition and core size. Unfortu- 
nately, these models require one to make assumptions 
about the degree of differentiation, among other things 
(|Baraffe et al.l 12008). A good measurement of k 2p , how- 
ever, reveals important independent structure informa- 
tion, which can break the degeneracies between bulk 
composition and the state of differentiation. Given such 
a wide range of potential k 2p values, even an imprecise 
measurement of k 2p will be extremely valuable for un- 
derstanding extra-solar planets. By measuring the k 2p 
values for extra-solar planets, we can also uncover con- 
straints on the density structure that are independent of 
the measurement of the planetary radius. This new in- 
formation may allow us to probe the unknown physics 
responsible for the currently unexplained radius anoma- 
lies. 

2.2.1. Induced External Gravity Field 

The internal structures of planets in our own solar sys- 
tem are most readily characterized by the zonal harmon- 
ics of the planet's gravity field, i.e. J 2 , J4, etc. It is 
these high-order harmonics that are directly measured by 
spacecraft flybys. To better understand the connection 
between the two, we can relate the k 2 formulation to J 2 
by writing out the expression for the induced potential at 
the surface of the planet in Equation[T]in terms of the def- 
inition of J 2 , yielding: k 2p V 2 app (R p ) = — J 2 —j^-P 2 {cos 8), 

where P 2 is the usual Legendre polynomi al and 9 is the 
planetary co-latitude (Murr ay fc Dermott 1999). We can 
use this equation to obtain expressions for the J 2 field 
induced by both rotation and tides (discussed in more 
detail below). The relation relies on dimensionless con- 
stants which compare the strength of the acceleration 
due to gravity with that of the rotational and tidal po- 
tentials: 
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v 2 R 3 
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where v p is the angular spin frequency of the planet. 
For the case where the spin axis and tidal bulge axis 
are perpendicular (i.e. zero obliquity), the relationship 
between J 2 and k 2 is, to first order: 



J-2 



(Jr 



(3) 



Note that q t is a function of the instantaneous orbital 
separation, r, and is thus constantly changing in an ec- 
centric orbit in response to the changing tidal potential. 
Hence J 2 for eccentric extra-solar planets is a complex 
function of time. This is why it is more sensible to ana- 
lyze the orbital precession in terms of k 2 , which is a fixed 
intrinsic property of the planet, rather than J 2 . 

As very hot Jupiters are expected to be synchronously 
locked (denoted by s) with small eccentricities, it can eas- 
ily be shown that g| w — 3q r , which simplifies equation[3] 
yielding: 
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(4) 



Using a moderate value of k 2p = 0.3, the J 2 of very hot 
Jupiters reaches as high as 5 xlO -3 , about half of the 
measured J 2 of Jupiter and Saturn. 
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2.3. Apsidal Precession 

The quadrupole field created by rotational and tidal 
potentials discussed above induces precession of the star- 
planet orbit. Both Jupiter and Saturn have rather signif- 
icant quadrupoles, dominated entirely by their sizeable 
rotational bulges resulting from rapid rotation periods 
of less than 10 hours. In contrast, very hot Jupiters are 
expected to be synchronously rotating, and thus their 
spin periods are longer by a factor of a few. Since the 
rotational bulge size goes as the square of the spin fre- 
quency, very hot Jupiters should have rotational bulges 
that are at least an order of magnitude smaller than 
Jupiter and Saturn, inducing only tiny quadrupole fields. 
These extra-solar planets are extremely close to their par- 
ent stars, however, with semi-major axes of only ~ 6 
stellar radii. Very hot Jupiters are thus expected to have 
large tidal bulges which are shown below to dominate the 
quadrupole field and resulting apsidal precession. 

2.3.1. Precession Induced by Tidal Bulges 

The orbital effect of tidal bulges is complicated by their 
continuously changing size. While tidal bulges always 
point directly 4 at the tide-raising object, their size is 
a function of orbital distance. Since the height of the 
tidal bulge depends on the actual separation between 
the objects, the second-order gravitational potential is 
time- varying in eccentric orbits. Accounting for this de- 
pendence (which cannot be capture d by us ing a fixed J 2 ) 
is critical, as illustrated by Sterne (1939a|). The domi- 
nant tidal perturbation to the external gravity field of 
the planet, evaluated at the position of the star, is a 
second-order potential: 



1 



V tid (r) = -k 2 GM*R b p r- 



(5) 



The apsidal precession due to the tidal b ulge, including 
the effect of bot h the st ar and the pla net is (|Sternelll939at 
lEggleton k, Kiseleva-Eggleton|[200l[) : 

Wtidal = C^tidal,* + k>tidal,p 
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where n is the mean motion and /2(e) is an eccentricity 
function: 



/ 2 (e) = (l- e 2 r 5 (l + ^e 2 + i e 4 ) 



1 



13 
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(7) 



Note that the factor of 15 does not appear for stationary 
rotational bulges, as detailed below, and comes through 
Lagrange's Planetary Equations from the higher depen- 
dence on radial separation (r~ 6 ) in the tidal potential. 
For this reason, tidal bulges are much more important in 
producing apsidal precession. 

Furthermore, the main factor of importance to extra- 
solar planets is the mass ratio, which comes in because 

4 We can ignore the lag due to dissipation, which has an angle 
of only Qp 1 < 10~ 5 for giant planets (Gold reich fc Soterl 119661 : 
[Murray & Dcrmott 1999T). 



the height of the tide is proportional to the mass of the 
tide-raising body. Consider the ratio of the planetary 
and stellar effects: 



Wtidal,p _ K 2p 
Wtidal.* &2* 
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For tidal bulges, the apsidal motion due to the planet 
clearly dominates over the contribution of the star. Even 
though the planet's radius is smaller than the star's by 
a factor of ten, the star is so much more massive than 
the planet that it raises a huge tidal bulge, which con- 
sequently alters the star-planet orbit. The benefit pro- 
vided by the inverse square of the small mass ratio is 
compounded by the order of magnitude increase in k 2 of 
the planet over the star. 

2.3.2. Precession Induced by Rotational Bulges 

The quadrupolar gravitational field due to the plan- 
etary rotational bulge, evaluated at the star's position 
is: 

V Iot (r) = h 2 v p l R 5 p r- 3 P 2 (cosa p ) (9) 

where a p is the planetary obliquity, the an gle between the 
orbit normal and the planetary spin axis. ISternd (|1939a|) 
assumes zero obliquity and calculates the secular effect of 
this perturbation on the osculating Keplerian elements. 
This final result, including the effect of both the star and 
the planet is 5 : 
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where g 2 (e) is another eccentricity function: 



g 2 (e) = (l-e 2 )- 



l + 2e 2 



3e 4 



(10) 



(11) 



Evaluating the importance of this effect requires an un- 
derstanding of the spin states of very hot Jupiters and 
their stars. The rotation and spin pole orientation of very 
hot Jupiter s should be tidally damped on timescales < 1 
MYr (e.g., iDobbs-Dixon et all 120041 : iFerraz-Mello et all 
2008). We therefore assume that all planets have rea ched 
the p suedosynchronous rotation rate derived by iHutl 
(1981). The rotation rate of the star is usually much 
slower since the tida l stellar spin-up time scale is much 
longer than -1 GYr (jFabrvckv et alj|2007l ). 

If both the star and the planet were spinning syn- 
chronously, the stellar and planetary rotational bulges 
would have comparable contributions to apsidal preces- 
sion. However, since the tidal bulge of the planet is 
a much more important effect, we find that even fast- 
spinning stars have a very weak contribution to apsidal 
precession. 

5 The f ull equ ation, including arbitrary obliq uities, is given 
IKopall $Wm, Equation V.3.1 8 (see also ISternd I1939al: 



Egglcton & Kisclcva-Egglcton 2001). Also recall that, unlike these 
authors, we use the symbol ki to represent the Love number which 
is twice the apsidal motion constant called k^ in eclipsing binary 
literature. 
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2.3.3. Total Apsidal Precession 

The other major contributor to the apsidal precession 
in extra-solar planetary systems is general relativity. The 
anomalous apsidal advance of Mercury's orbit due to its 
motion near the massive Sun was one of the first con- 
firmations of general relativity. This same apsidal ad- 
vance is prevalent in very hot Jupiter systems and has 
been shown to be possibly detectable through lon g-term 
transit timing (Miralda-Escudc 2002; Hcvl & Gladman 
120071 : iPal fc Kocsid 120081 : IJoTdan fc Bakod I2008D . The 
relativistic advance is given (to lowest order) by: 



^GR 



3GM*n 

ac 2 (l — e 2 ) 



(12) 



One additional eff ect for non-synchronou s planets is 
due to thermal tides (| Arras fc Socratesll2009D . which cre- 
ate a bulge on the planet due to temperature-dependent 
expansion of an unevenly-radiated upper atmosphere. 
The thermal tidal bulge is very small in mass and is not 
expected to provide a significant contribution to apsidal 
precession (P. Arras, pers. comm.) and is thus neglected. 

Since we are considering only the lowest-order effects, 
all the apsidal precession rates (rotational/tidal for the 
star/planet and general relativity) simply add to give the 
total apsidal precession (roughly in order of importance 
for very hot Jupiters): 



<^tot — Wtid,p + WGR + W rot! p + UJ m t,* + Wtid, 



(13) 



We are ignoring the small cross-terms (geodetic preces- 
sion, quadrupole-quadrupole coupling, Lense-Thirring 
effect, nutation, etc.) for the purposes of this paper as 
higher-order corrections. 

Calculating each of these contributions to the preces- 
sion shows that for very hot Jupiters, the dominant term 
in the total apsidal precession is due to the planetary tidal 
bulge. For the known transiting planets, the fraction of 
apsidal precession due to the planet is calculated and il- 
lustrated in Figure 1. The precession due to the interiors 
of very hot Jupiters towers over the other effects. Gen- 
eral relativity, the next largest effect is ~10 times slower 
than the precession caused by the planetary tidal bulge. 

The apsidal precession rate of very hot Jupiters due 
solely to the interior structure of the planet is: 



3.26 x 10~ 1U rad/sec x 
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(14) 



which explains why low density very close-in Jupiters are 
the prime targets for measuring apsidal precession. For 
these planets, the precession rate can reach a few degrees 
per year. 
The precession due to the planet has generally been ne- 
lected in extra-s olar p l anet transit timing work to date 
Miralda-Escudi 120021 : iHevl fc GladmarJ I2007T ). which 
has considered stellar oblateness or general relativity to 
be the domi nant effects (in the abs ence of other plan- 
ets) though Uordan fc Bakosl ((2008) have also pointed 
out that Wtidai.p can be an important source of apsidal 
precession. We find that the planetary quadrupole is 
usually 1-2 orders of magnitude more important than ef- 
fects previously considered for single very hot Jupiters. 



l 

0.9 

0.8 

§ 0.7 

| 0.6 

| 0.5 

1 0.4 

| 0.3 

0.2 

0.1 





♦ + 



♦ ♦ 

♦♦ + 



2 3 4 

Period (days) 



Fig. 1. — Fraction of Apsidal Precession Due to the Plan- 
etary Quadrupole. The points show the planetary fraction of 
the total apsidal precession calculated for the known transiting 
extra-solar planets with properties taken from J. Schneider's Extra- 
Solar Planet Encyclopedia (http://www.exoplanet.eu), assuming 
the planet has a typical Love number of ki v = 0.3 (e.g. Saturn- 
like). The apsidal precession induced by the tidal and rotational 
bulges of the planet overcome precession due to general relativity 
and the star, especially for short period planets. The "error bars" 
show the range of planetary contributions for a 5% variation in 
stellar masses (and hence tligj) and the comparatively smaller ef- 
fect of varying the stellar Love number and rotation rate over all 
reasonable values. The five cases where the planetary contribu- 
tion to apsidal precession is most important (boxed) also have the 
shortest precession periods: WASP-12b, CoRoT-lb, OGLE-TR- 
56b, WASP-4b, and TrES-3b would fully precess in about 18, 71, 
116, 120, and 171 years, respectively. The planet in the lower left 
is CoRoT-7, a super-Earth planet whose planetary contribution to 
precession is small because of its small radius. Transiting planets 
with periods longer than 6 days all had planetary contributions less 
than 0.15. In all cases, the dominant signal in apsidal precession 
of very hot Jupiters is fep, which is determined by their inter- 
nal density distribution and is a powerful probe into their interior 
structure. 



Hence, measuring apsidal precession essentially gives 
^tid.p which is directly proportional to fep, implying that 
transit light curve variations due to apsidal precession 
can directly probe the interiors of extra-solar planets. 

2.4. Modification of the Mean Motion 

Non-Keplerian potentials also modify the mean- 
motion, n, and cause a small deviation from Kepler's 
Third Law. Including the effects described above, the 
non-Keplerian mean motion, n', is (dropping second- 
order corrections): 



= n 1 



3GAf» 
2ac 2 
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where e is defined as 
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and n 2 _ Gm 



The general relativistic correction to 
the mean motion is from lSoffell Jl989). (Throughout this 
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paper, except where noted, the difference between n' and 
n is ignored as a higher-order correction.) 

As with apsidal precession, the planetary quadrupole 
is more important than the stellar quadrupole by about 
2 orders of magnitude. At the largest, t he cor r ection 
to the mean motion is a few times 10~ 5 . Iloriol (|2006l ) 
used the fact that quadrupole moments cause deviations 
to Kepler's Third Law to attempt to derive the J2 of 
the star HD 209458 (the quadrupole of the planet was 
incorrectly ign ored) . 

However, as Iloriol ()2006l ) found, this method is only 
feasible if you know the masses and semi-major axes 
of the orbit a priori or independently from Kepler's 
Law. Since the error in stellar masses (from radial ve- 
locities and evolut ionary codes) is usually 3-10 % (e.g. 
iTorres et al.ll2008ft . the propagated error on k,2 V would 
be a few times greater than the highest /c2 P expected, 
making this method impractical. It has been proposed 
that the stellar mass and semi-major axis can be pre- 
cisely and independentl y measured via the light-travel 
time effect described by lLoebl J2005). in practice, how- 
ever, the light-travel time effect is highly degenerate with 
the unknown transit epoch and/or the orbital eccentric- 
ity. We find that a precise independent measurement of 
M* from light-travel time is impractical even with the 
excellent photometry of Kepler. 5 

2.5. Expectations for Planetary Eccentricities 

Thus far, we have quantified how planetary interiors 
affect the orbit through precession. The photometric ob- 
servability of this apsidal precession is highly dependent 
on the current orbital eccentricity (e). Small eccentrici- 
ties are the largest limitation to using transit light curves 
to probe extra-solar planet interiors. Indeed, if eccen- 
tricities are very low, measuring apsidal precession from 
transit light curves may not be possible for any of the 
Kepler planets. 

Nearly all hot Jupiters have eccentricities consistent 
with zero, though the radial velocity technique has dif- 
ficulty putting 3-g upper limits on eccentricities smaller 
than 0.05 (|Laughlin et al. 2005). So far, the strongest 
constraints are placed by comparing the deviation of the 
secondary transit ti me from half the orbital pe riod, which 
are related by (e.g.. lCharbonneau et al.ll2005l ): 
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Similarly, by measuring the primary and secondary 
transit durations (0/ and On), an additional constraint 
can be placed on e sin w. The equati on commonly quoted 
in the extra-solar p l anet literature (iKallrath et aLlll999l : 
ICharbonneaul 120031 : 1 Winn et"ai1 l2006ft has a s ign error; 
the correct equation is derived byjKopal (1959), p. 391 : 
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where a = — * ^ . The accuracy of this measurement 

a^/l — e 2 

is typically smaller than for ecosw, but we include this 

6 We do note that detailed observations of multiple-planet sys- 
tems can yield mass estimates of each of the bodies independently. 
Kepler asteroscismology can also provid e independent informa tion 
about stellar mass and other properties (Kjeldscn et al. 2008). 



equation to note that there is information about both the 
eccentricity an d its orientation i n the full transit light 
curve (see also lBakos et al. 2009). 

Combining secondary transit timing information with 
radial velocity and Rossiter-McL aughlin measurements 
to help constrain u>, I Winn et"ail (|2005ft found the best- 
fit eccentric i ty fo r HD 209458 was ~0.015. Though 
IWinn et al.l (|2005f ) argue that the actual ecce ntricity is 
proba bly less th an 0.01, it is not n ecessarily (Mardling 
l2007h . Recently. IJoshi et alJ (J2008ft revealed WASP-14b, 
a young massive hot Jupiter with an eccentricity of 0.1; 
WASP-lQb and WASP- 12 b also appear to be eccentric 
(|Christian et all 120081 : iHebb et alj |2009). though these 
eccentricities may be spurious or overestimated. 

The most accurate eccentri city constraint is a detec- 
tion by iKnutson et al.l (|2007al ) for the very hot Jupiter 
HD189733b. They observed continuously and at high 
cadence (0.4 seconds) with the Spitzer space telescope 
and measured a secondary timing offset corresponding 
to e cos uj — 0.001 ± 0.0002, a 5-cr result that they could 
not explain by any other means. (Preliminary analy sis of 
additional data for this planet bv lAgol et al.l (J2009) indi- 
cates ecosu = 0.0002±0.0001.) The constraint one sin w 
is much weaker. A non-zero eccentricity of e ~ 0.003 for 
hot Jupiters is therefore consistent with every measure- 
ment available in the literature, though the actual values 
of eccentricities at the 10 ~ 3 level are essentially uncon- 
strained. 

In the absence of excitation, the current eccentricities 
of these planets depend on the initial eccentricity and the 
rate of eccentricit y decay. Extrapolating f rom planets in 
our solar system ()Goldreich fe Soterlll966f) implies short 
circularization timescales of ~ 10 MYr, though recent 
studies have shown that using a fixed eccentricity damp- 
ing timescale is an in appropriate simplification of the full 
tidal evolution (e.g. iJackson et "all 120081 : iLevrard et al.1 
I2009t iRodriguez fe Ferraz-Mellol 12009ft . Even an analy- 
sis using the full tidal evolution equations cannot give a 
compelling case for the present-day eccentricities of these 
planets, since there are essentially no direct constraints 
on the tidal dissipation parameter for the planet, Q p . 
Various estimates show that Q p for exoplanets is not 
known and may be quite large (e.g., iMatsumura et all 
2008), implying that non-zero eccentricities are not im- 
possible. Even so, we stress that the best candidates for 
observing apsidal precession are also those planets that 
have the fastest eccentricity damping, since the damping 
timescale and apsidal precession rates are both propor- 
tional to jf- (jr) ■ Hence, those planets which have 

the fastest precession rates will also have the lowest ec- 
centricities. The first step in determining if this trade-off 
allows for apsidal precession to be measured by Kepler 
data is to apply the techniques described in this paper 
to the data themselves. Furthermore, with the discov- 
ery and long-term characterization of more planets using 
ground and space-based observations, the detectability 
of apsidal precession will increase dramatically. 

We should note that there are several mechanisms 
that can excite eccentricities and compete with or over- 
whelm tidal dissipation. The most prevalent is as- 
sumed to be eccentricity pumping by an additional 
companion dPeale et al.1 119791 : iBodenheimer et aTll2001t 
I Adams fc Laughlm1l2006f ). Even very small (Earth-mass 
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or less) companions in certai n orbits can pro vide signif- 
icant eccentricity excitation (jMardlind [20071 . (In this 
case, however, our single-planet method for estimating 
fep would need to be modified considerably.) Tidal dis- 
sipation in rapidly rotating stars tends to increase the ec- 
centricit y, potentially prolonging circularization in some 
systems (jFerraz-Mello et aLll2008t ). Very distant inclined 
companions (e.g. a planet orbiting a star in a mis- 
aligned binary star system) can induce Kozai oscillations 
that i mpart very large eccent r icities on secular timescales 

e.g., iFabrvckv fc Tremainel l2007ft . lArras fc Socrates! 

2009) proposed that thermal tides can significantly af- 
fect the orbital and rotational properties of extra-solar 
planets , though th eir c onclusions app ear to be overesti- 
mated (iGoodmanll2609HGu fc Qgilvieil2009T) . Finally, re- 
cent (not necessarily primordial) dynamical instabilities 
in the planetary system can also be responsible for gener- 
ating eccentricity which simply hasn't damped away yet 
dFord et all 120051: I Gomes et al.l [2005 IChatteriee et all 
l2007tlThommes et al.ll2008|) . We. therefore, continue our 
analysis under the possibility that some very hot Jupiters 
may have non-zero eccentricities. 

3. TRANSIT LIGHT CURVES OF APSIDAL PRECESSION 

Previous studies of transit light curve variability due 
to non-Keplerian perturbations have focused almost ex- 
clusively on transit timing. In contrast, we model the 
full photometric light curve in order to estimate the de- 
tectability of k2 P . This will automatically include the 
effect of changing transit durations, which are very use- 
ful for detecting apsid al precession ()Pal fc Kocsisll200"8l ; 
Uordan fc Bakosl[2008D . In addition, using full photom- 
etry can provide a more direct and realistic estimate 
of the detectability of h^-p. Of course, the drawback is 
additional computational cost, though we found this to 
be manageable, requiring less than 20 seconds to gener- 
ate the ~ 2 million photometric measurements expected 
from Kepler's 1-minute cadence over 3.5 years. 

3.1. Our Transit Light Curve Model 

Determining the photometric light curve of a transit- 
ing system requires knowing the relative positions of the 
star and the planet at all times. These can be calcu- 
lated by describing the motion of the planet with time- 
varying osculating orbital elements. When describing 
the motion of the planet using instantaneous orbital el- 
ements, it is usually customary to ignore the periodic 
terms by averaging, as in Sterne (1939a), and calculate 
only the secular terms. These small periodic terms de- 
scribe how the orbital elements change within a single 
orbit as a function of the true anomaly, /, due to the 
non-Keplerian potential. In precessing systems, the value 
of the true anomaly at central transit, f tr = 90° — uj tr , 
changes subtly from one transit to the next, inducing 
slow variations in the osculating orbital elements at tran- 
sit. Therefore, we include in our model the dominant 
periodic changes in orbital elements as a function of or- 
bital phase, using M tr ~ ftr as an appropriate approx- 
imation for low eccentricities. Using a direct integra- 
tion (described in Section I4TTJ) . we verified that ignoring 
these periodic variations can cause non-negligible sys- 
tematic errors in determining transit times. The peri- 
odic changes are derived from the same disturbi ng po- 
tentials used above. We follow the method of iKozail 



(1959) for calculating osculating elements from mean el- 
ements, and assume zero obliquity. The correction is 
similar to the correction to the mean motion, which is 
also applied in our model. The correction to the semi- 
major axis, eccentricity, longitude of periapse, and mean 
anomaly are a osc = a moa n + pr% e cos M « 2aeecosM, 

e 0SC = e mcan + e(l ~ COSM), W osc = W moan + fshlM, 

and M osc = M moan — - sin M where e is defined in Equa- 
tion [THJ General relativistic period i c corr ections are also 
added; these are taken from Soffel (1989), page 92 (with 
a = 0, /3 = 7 = 1). Using our direct integrator (de- 
scribed below), we verified that these corrections repro- 
duced the actual orbit to sufficient accuracy for this anal- 
ysis as long as e > e ~ 10~ 5 . Other corrections are 
higher order in small parameters and are ignored. 

Our model uses these corrected elements to generate 
astrocentric Cartesian coordinates for a specific system 
inclination and, for completenes s, also includes the effect 
of light-travel time (lLoebll2005h though we con cur with 
Uordan fc Bakoi (|2008l ) and lPal fc Kocsisl (|2008D that the 
light-travel time change due to u is unimportant. The 
positions are then translated to photometric light curves 
using the quad r atic li mb-darkening code 7 described in 
iMandel fc Ago! (|2002t ). Kepler data will have enough 
signa l-to-noise to justify usi ng non-linear limb darkening 
laws ([Knutson et al.ll2007rj) . but we do not expect that 
this simplification will significantly alter our conclusions. 

In addition, we include the photometry of the sec- 
ondar y eclipse. As suggested bv lLopez-Morales fc Seagerl 
(|2007| ). very hot Jupiters can reach temperatures ex- 
ceeding 2000 K, where their blackbody emission at op- 
tical wavelengths is detectable by Kepler. This thermal 
emission is added to the reflected light of the planet, 
which appears to be small based on the low upper limit 
of the albedo of HP 209 458b and TrES-3 measured 
by iRowe et al.l (|2007l ) and IWinn et all (|2008f l. respec- 
tively. We fin d that in Kepler 's observing bandpass 
of 430-890 nm (jKoch et al.ll2006( ). thermal emission of 
very hot Jupiters can dominate over the weak reflected 
light. We estimate the depth of the secondary eclipse 
(d SGC ) in our simulated Kepler data by assuming that 
1% of the light is reflected and the other 99% absorbed 
and reemitted as processed thermal blackbody emission 
from the entire planetary surface (day and night sides). 
To be conservative and to account for unmodeled non- 
blackbody effect s, we divide the resulting planet/star 
flux ratio by 2 (|Hood et al.l l2008f) ; the resulting depth 
of around 2 x 1 0~ 4 is consistent with the lower values of 
iBurrows et al.l (|2008f ). the tentative measurement of th e 
thermal emission from CoRoT-2b (jAlonso et al.l 12009) . 
and the detecti on of secondary eclipse em ission from 
OGLE-TR-56b (jSing fc Lopez-Morales! l2009h . We note 
that the best candidates for detecting &2 P are those with 
small semi-major axes and large radii; these same plan- 
ets have relatively large d scc values (Table 1). Secondary 
eclipses are very useful for determining e and u. We will 
also find that they can be important for observing apsidal 
precession. 

Our model generates accurate photometry for an extra- 
solar planet undergoing apsidal precession. Several other 
small photometric effects have been discussed in the lit- 
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erature, which we do not include. Most of these effects 
are periodic (e.g. the reflected light curve) and therefore 
will not affect the long-term trend of precession. Care 
will need to be taken to ensure that slow changes due to 
parallax and proper motion, which should be qui te small 
for re latively dista nt stars observed by Kepler (jRafikovl 

f)08; Scharf 2007) or changes in the stellar photosphere 
joeb 2008) are not significant. Non-Gaussian astro- 
physical noise of the star and other systematic noise 
should degrade the accuracy with which k 2p can be mea- 
sured compared to our ideal photometry. The long-term 
varia bility of the star c an be interpolated away or mod- 
eled (|Lanza et a l. 2009), though it is not clear how short- 
term variability will affect transit light curves at Kepler's 
level of precision. On the other hand, complimentary 
observations (e.g., warm Spitzer, HST, radial velocities, 
JWST, etc.) should only enhance our understanding of 
the systems studied. 

3.2. Accuracy of k 2p measurement 

With an accurate photometric model of apsidal pre- 
cession, one could estimate the measurement accuracy of 
k 2p from Kepler data by carrying out a full Monte Carlo 
study of the inversion problem, going from realistic syn- 
thetic photometric data sets to a determination of all 
system parameters. In this work, instead, we carry out 
a much simpler calculation which cannot provide strict 
one-sigma error estimates like the Monte Carlo analy- 
sis, but does give an indication of how well k 2p can be 
resolved given a large dataset. 

We obtain this accuracy estimate by comparing a re- 
alistic precessing photometric model with k 2p ^ to a 
base model with k 2p = 0. The base model is still under- 
going very slow apsidal precession, induced by general 
relativity and k 2 * • We calculate the effect of a non-zero 
k 2p value by subtracting the precessing model from the 
base model. (See Figures [2] and |U) Then, by calculating 
the root-sum-square of the residual signal and comparing 
it to the photometric error on a single data point, we ob- 
tain a numerical measure of the relative signal induced 
by k2 P . The "signal-to-noise" ratio for the data set is 
therefore given by: 



1 
N 



VEM^W 



(19) 



where yi and y® are the photometry model values for the 
k 2p test model and the base model, respectively, and a is 
the photometric error. We use a — 1000 parts per million 
(ppm) flux per 1-minute integration, corresponding to 
the expected nois e of Kepler on a faint V — 14 star 
(jKoch et alJl200a ). Of the 30 planets with periods less 
than 3 days, 16 are expected to be brighter than V ~ 14 
(T. Beatty, pers. comm.) and we can reasonably expect 
some fraction of these to have orbits comparable to the 
planets modeled here. 

Since our residual signal changes as a function of time, 
this is not a true signal-to-noise calculation; the distri- 
bution of values in time matters for a proper interpreta- 
tion, but any distribution would yield the same effective 
jt , and thus this construction is not capturing all of the 
details. Even so, it does provide a useful and reasonable 
rough estimate for detectability. In order to identify the 
resolution on the k 2p measurement, we search for the 
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Fig. 2. — Photometric Difference Signal from k2 p . As 
described in the text, we use the difference between two theoretical 
light curves in the transit photometry to assess the observability of 
apsidal precession by Kepler. For WASP-4b at u = 0°, e = 0.003, 
and a central impact parameter, the difference between a model 
with A)2p = and ki v = 0.146 would yield an effective "signal- 
to-noise" of 1 on a moderately bright star (V = 14). Shown is 
this difference signal; the root sum of squares of the signal is equal 
to 1000 ppm, the expected photometric accuracy of Kepler for a 
1 minute observation (Koch et al. 2006). The trends seen in the 
figure are illustrated in Figure \3\ by considering excepts of single 
primary transits from the regions labeled 1-5. 



value of k 2p which yields a signal-to- noise of j? = 1. 
This is reasonable since it represents the threshold value 
of k 2p , below which planetary induced precession cannot 
be distinguished in the data with the given errors. The 
threshold k 2p value can also be loosely thought of as an 
estimate of the 1-cr expected errors. 

This is a realistic estimate only insofar as the resid- 
ual signal {yi — y®) is due only to k 2p and cannot be 
absorbed by any other parameters. Hence we seek to 
choose other parameters so as to minimize the residu- 
als without changing k 2p . For most system parameters, 
this is accomplished by referencing the time to the cen- 
ter of the data set, and thus the difference between the 
signals grows similarly forward and backward in time as 
seen in Figures [2] - [5] The transit shapes in both models 
are equivalent at the center of the dataset as would be 
expected in an analysis of actual data. 

Additionally, a major effect from changing the preces- 
sion period is to alter the observed average period. When 
analyzing actual data, this would just be absorbed into a 
small adjustment to the (unknown) stellar mass, thereby 
adjusting the period to absorb much of the k 2p signal. It 
is therefore important to correct for the average period 
change to avoid significantly overestimating the signal 
due to k 2p . Additionally, there is a similar, though less 
severe, effect for the epoch of the first transit, which is 
also adjusted to best absorb signal. This is achieved by 
using an analytic expression for the transit times (see 
Equation |2"21 below) which match the transit times of the 
photometric model to very high accuracy. By fitting a 
line to these times, we can determine the average period 
and epoch that absorb the degenerate portions of the k 2p 
signal, leaving behind the residual due only to k 2p . We 
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Fig. 3. — Excerpts of Photometric Difference Signal. Ex- 
amining excerpts of the residual signal shown fully in Figure [2] the 
effects of transit timing and "transit shaping" can both be seen. 
The five excerpts are offset for clarity. Transit timing has an asym- 
metric signal (dotted lines), obtained when subtracting two transit 
curves slightly offset in time. Transit shaping, which is mostly due 
to changing transit duration, creates a symmetric signal (dashed 
lines). The total difference signal (solid lines) is dominated by the 
effect of transit shaping, which has ~30 times more signal than 
transit timing alone. (See explanation in text.) Both effects are 
maximized at the beginning (1) and end (5), as expected for a 
signal that increases with longer baseline. The maximal signal oc- 
curs during ingress and egress, when the light curve changes the 
fastest. The transit shapes are equivalent at the center (3) by con- 
struction. The transit timing anomaly of precession is quadratic, 
which, when fitted with a best-fit straight line corresponding to 
a non-precessing signal, yields two intersections when transit tim- 
ing is minimized (2,4). The transit timing offset at the beginning 
and end is only 0.085 seconds, while the center is offset by -0.042 
seconds. 



have not explicitly accounted for degeneracies between 
the signal from k^p and the other parameters, like the 
radius, limb darkening, and system inclination, but since 
A;2p induces a time varying signal while these other pa- 
rameters are generally constant, there is little expected 
signal absorption from these parameters. 

The only major drawback of this approach is that it 
does not allow the eccentricity state of the system to 
change. With real data, the eccentricity and preces- 
sion phase are not known in advance, and thus must 
be found by inversion. As detailed in Section l2~5l eccen- 
tricity and orbital orientation are primarily constrained 
by comparing primary and secondary transit pairs, and 
thus proper inversion is greatly aided by accurate ob- 
servations in wavelengths more favorable to secondary 
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Fig. 4. — Photometric Difference Signal from k2 p . Similar 
to Figure [2] but for lo = 90° . This figure is dominated by the 
photometric difference between secondary transits slightly offset 
in time. At u = 90° the changes in the primary transits due 
to precession are small, except far away from the central time. 
At this orientation, the primary-secondary timing offset (Equation 
1 1 76 is maximized. This "secondary transit timing" signal is weaker 
than the signal from primary transit as the secondary transit depth 
is much shallower. Therefore, an unreasonably high ki v of 0.925 
is required to detect the apsidal precession. Excerpts of single 
secondary transits taken from regions labeled 1-5 are shown in 
Figure \E\ 

transit observat ions, obtained by Spitzer, HST, or from 
the ground (e.g.lKnutson et al.|[2007ilSwain etaill2008l : 
iGillon et al.ll2008l) . We also find that binned and folded 
Kepler data has comparable sensitivity to a single Spitzer 
observation for characterizing the secondary eclipses of 
very hot Jupitcrs. In any case, our assessment of the 
threshold k^ v assumes that the eccentricity of the system 
is very well known, which will likely require additional 
supporting observations. 

3.3. Comparison to Expected Signal 

The residual light curves calculated for each planet, 
Figures [2] - El match the theoretic al expectations 

of the apsidal preces si on s ig nal (IMiralda-Escudel 

2001 iHevl fc Gladmanl [20071 : iPal fc Kocsisl 120081 : 



Jordan fe Bakosll2008f) . To interpret the results of our 



analysis, it will be useful to briefly review the major 
components of the apsidal precession signal: changes in 
the times of primary transits, changes in the shape of 
primary transi ts, and changes i n the p rimary-secondary 
offset times (IMiralda-Escudel l2002t iHeylfc Gladmanl 
[2rM[PlTfc Kocsisll200a Uordan fc Bakosll2008ft . " 

The primary transit times, T/v, due to apsidal preces- 
sion are well described by a sinusoid for very low eccen- 
tricities (e <C 0.1): 



Ti 



N 



To 



'-(cOSUJtr.N — COSU)tr,o)] (20) 



where Tq is the epoch of the first transit, 
u> tr ,N — w(Tjv — T ) + u> tr o is the argument of peri- 
apse for the N th transit, and P b s is the observed period 
between successive transits, which deviates from the 
actual orbital period since the o rbit ha s proc essed a 
small amount between transits (|Battenl Il973f) . For 
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Fig. 5. — Excerpts of Photometric Difference Signal. Sim- 
ilar to Figure l3l but for w = 90°. Single secondary transit differ- 
ences are excised from the full difference signal shown in Figure 
3] The shape of the curves is due to the subtraction of two sec- 
ondary transits slightly offset in time. Since the secondary tran- 
sits are complete occultations, they are flat-bottomed and lack the 
additional structure due to limb-darkening seen in Figure [3] By 
construction, the offset grows in time away from the center (3) of 
the signal and attains a maximum at the beginning (1) and end 
(5). Curves 2 and 4 are shown for comparison to Figure [3] 

small eccentricities, the amplitude of the transit timing 
variations due to ki v is: 



eR 
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119 sec x 



V0.003/ U.025 AU/ 



3/2 / M* 



\M C . 







-1/2 



( 21 ) 

Given that individual transit times can be measured with 
accuracies of only a few seconds, even tiny eccentricities 
e < 10 -5 can induce detectable transit timing variations 
on precessional timescales (~ w" 1 )- 

For our analysis, we extended Equation [20] to fifth or- 
der in eccentricity allowing accurate determination of 
transit times for eccentricities up to of order 0.1. We 
also require a correction for the effect of a non-central 
impact parameter (i < 90°, e > 0). For an inclined 
eccentric orbit, the apparent path of the planet across 
the stellar disk is curved. At orientations where the line 
of sight is not along the major axis of the ellipse, the 
curved path is also asymmetric. Therefore, the times of 
photometric minima, T/v, do not correspond exactly to 
the times of conjunction (when the planet crosses the 
y — z plane and f tr = 90° — u^ A We follo w the cor- 
rection from Equation VI. 9-21 of Kopal (1959), who find 



that at photometric minimum, f tr = 90° — u>' tr , where 
u>' tr = u> tr + e cos u> tr cot 2 (i) (I — e sin u> tr esc 2 (i)); in this 
corrective term, it is only required to keep terms up to 
second order in eccentricity. Assuming that i and u> are 
constant, it can be shown that 

T N = T + NP ohs 

, -*obs 



e{cosu>' t N - cosu' t ) 



3 

-e (sm2u)L N - sin2u/L ) 



-e 3 (cos3Wf 



; ,, v -cos3wL ) 



^[j^(sin2uj' trN - sin2o4,o) 
5 



cos 3lu', 



5 3w tr _ ) 



_3_ 

80 



(cos5a4_jv-cos5a4 )0 )J 



(22) 



This transcendental equation is solved iteratively for 
(Tjy — To) to obtain the transit times and has been tested 
thoroughly against the empirical determination of tran- 
sit times calculated by our light curve model described 
above. 

The expected apsidal precession periods (including 
small contributions from GR and the star) for WASP- 
12b, CoRoT-lb, OGLE-TR-56b, WASP-4b, and TrES- 
3b are around 18, 71, 116, 120, and 171 years, respec- 
tively. In other words, they have precession rates induced 
by the planetary tidal bulge of a few degrees per year, 
compared to a few degrees p er century as the fastes t gen- 
eral relativistic precession (|Jordan fe Bakosl 12008) . We 
caution that if — E for WASP- 12b is overestimated due 
to imprecise data (e.g. IWinn et"aT]|2007D . then the pre- 
cession period would increase accordingly. 

Even with such fast precession rates, the duration of 
observations will generally be much shorter than the pre- 
cession period. In addition, as discussed above, the 
linear timing anomalies will be absorbed into the ef- 
fective per iod as a small change i n the unknown stel- 
lar mass (iHevl fc Gladmanl [20071 iPal fc Kocsisl [20081 : 
iJordan fc Bako s 2008). Therefore, detection of apsidal 
precession from primary transit times alone will require 
a significant detection of the curvature over a small por- 
tion of a long-period sinusoid. Since the curvature in 
Equation [20] is maximal at u> rs 0, 180°, these orienta- 
tions have the best primary transit timing signal. Even 
at these orientations, detecting ki v from primary transit 
times alone is difficult, since it can be shown that the 
signal strength is proportional to eoj 2 , du e to the need to 
detect curvature (jHevl fc Gladmanl I2007T I. 

When the observational baseline is much shorter than 
the decades-long precession period, utilizing the chang- 
ing shape of the transits can signi ficantly improve de- 
tectability of apsidal precession (|Pal fc Kocsisl 120081 : 
IJordan fc Bakos 2001). Transit shapes are primarily de- 
termined by the orbital speed at transit ftr and impact 
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parameter 6, both of which depend on the precession 
phase uj tr . For small eccentricities, the orbital angular 
speed at transit is given simply by ft r — n(l + 2e cosuj tr )- 
Changes in the impact parameter are somewhat more 
subtle, since b is given by r tr cos i/R* , where r tr ~ 
a(l — e 2 )/(l + es'mujtr) is the star-planet separation. 
Hence, the apparent impact parameter of the planet can 
change for non-central transits, even when the orbital 
plane remains fixed. The evolving transit shape of pro- 
cessing orbits is determined by variations in both orbital 
speed and impact parameter. Simplifying the effect of 
transit shape by considering on ly the vari ations in tran- 
sit duration as a f unctio n of u>tr , iPal fc Kocsisl (J2008f ) and 
Uordan fe Bakosl (J20 08) find that these two effects are of 
comparable magnitude. These authors also show analyt- 
ically that the two effects exactly cancel when b = l/\/2- 
At this impact parameter, the transit duration stays con- 
stant throughout apsidal precession. The full photomet- 
ric transit shape, however, still changes detectably in a 
precessing orbit, though the magnitude of signal is re- 
duced (Figure [7|). 

The expected effect of changing transit shapes is fully 
consistent with the photometric difference signals calcu- 
lated by our model (Figures [2] and [3]). Indeed, our model 
shows that transit shaping dominates the signal by a fac- 
tor of >30 (Figure [3]). We can also see that changes 
in the transit shape are maximized at orientations near 
w~0, 180° (as expected from Equation [T5J), 

For small eccentricities, the transit shaping signal 



A 



oc eio oc eh 



2p- 



Therefore, 



strength is given by 

when transit shaping dominates the observable signal, 
we should find that searching for the threshold ki v value 
that yields ■£ = 1 results in a power law relationship 
between threshold ki v and e, such that /c2 P oc e . By 
solving for threshold &2 P for eccentricities from 0.001 to 
0.1, we find, as expected, that threshold ki v very closely 
follows a power law in eccentricity with a slope of -1 for 
all planets. This power law relationship can be written 
as ek2 P = C, where C is a constant calculated from our 
model that depends on the planetary, orbital, and stellar 
parameters of the system. 

At ijj w 90, 270°, transit timing and transit shaping ef- 
fects are much weaker and are rather ineffective at con- 
straining apsidal precession. At these orientations (when 
the Earth's line of sight is nearly aligned with the major 
axis of the orbit), another photometric signal emerges: 
variations in the difference between the times of primary 
and secondary transits. The changing orientation of the 
orbital ellipse causes a variation in the offset between 
primary and secondary transit times following Equation 
HZlabove (|Hevl fc Gladm an 2007t Uordan fc Bakosl2008l ). 
These authors show that the strength of this signal is 
also proportional to ecu and we find that the variation in 
threshold k2 P then also follows &2 P oc e _1 . 

The photometric difference signal at u> — 90° is shown 
in Figures |4] and El Using the method described in Sec- 
tion 13.21 to remove degeneracies almost eliminates the 
primary transit signal entirely, as expected, and the sec- 
ondary transit offset becomes the more powerful signal. 
For WASP-12b, with an expected Kepler secondary tran- 
sit depth of — 1830 ppm, the threshold ki v is actually 
lower at to = 90° (Figure [6]). For the other planets, the 
secondaries are not as important. 



Our estimates of threshold ki v at u> = 90° are based on 
the unknown secondary transit depth (d scc ) in the Ke- 
pler bandpass (though our estimates of d scc are consis- 
tent with all the measurements in the literature to date) . 
Furthermore, we find that 4> oc d scc , so that deeper sec- 
ondary transits improve the accuracy with which &2 P can 
be measured. It is important to note that combining Ke- 
pler primary transit times with precise secondary tran- 
sit times measured in the near-infrared (e.g. by warm 
Spitzer, HST, or JWST) is a very powerful way to con- 
strain apsidal precession (|Hevl fc Gladmanll2007t ) for any 
orientation. Even a few high-precision secondary eclipse 
observations are enough to lower the value of threshold 
/c2 P from our predictions, especially when 10 ~ 90, 270°. 

By construction, threshold k2 P values vary linearly 
with the assumed photometric error a = 0.001 x 
10°- 2 '^~ 14 ). In addition, re-performing our analysis us- 
ing a 6-year long Kepler mission improved threshold k2 P 
values by a common factor of —2.2. 

3.4. Results for Specific Planets 

Using the method described above, we have determined 
the threshold k2 P for the most favorable known transit- 
ing planets as analogs for the very hot Jupiters to be 
discovered by Kepler. The threshold fc 2p for each planet 
was computed at a range of eccentricities from 0.001 to 
0.1 and for u> = 0° and uj = 90°. Using the relation- 
ship discussed above (&2 P oc e _1 ) we interpolated (and 
sometimes extrapolated) our calculations to determine 
the eccentricity required to reach threshold /c2 P values 
of 0.3, 0.1, and 0.01. These results are summarized in 
Figure [6] and Table 1 . 

WASP- 12b is the best candidate for observing apsi- 
dal precession. With an eccentricity of e — 0.00026 and 
fc2 P =0.3, the apsidal precession would have an effective 
signal-to-noise of — 1 for all of Kepler data. If e is —0.001, 
then /c2 P can be well characterized and not just detected. 
As the difference in k,2 between Jupiter and Saturn of 
— 0.15 is primarily due to the presence of a massive core, 
a resolution in fc2 P of 0.1 is enough to detect whether or 
not the planet has a core, at the — 1-sigma level. 

Although WASP- 12b does not lie in the Kepler field, 
it clearly stands out as an excellent candidate for observ- 
ing apsid al precession . Though the putative eccentricity 
of 0.049 (lHebb et al.ll2009fl is probably an overestimate 
(jLaughlin et al.ll2~005h . if it were real, it would cause si- 
nusoidal transit timing deviations with an amplitude of 
—25 minutes (using Equation |2"0)) and a period of — 18 
years. Such a large deviation would be readily observed 
from the ground in either transit times or transit shapes. 
If apsidal precession is not observed, tight upper limits 
on the eccentricity can be established. 

Analogs to the very hot Jupiters WASP-4b, TrES- 
3b, CoRoT-lb, and OGLE-TR-56b are good candidates 
for observing apsidal precession if the eccentricities are 
above -0.003. (Note that CoRoT-lb has o nly -30 days 
of ob servations from the CoRoT satellite (|Barge et al.l 
l2008f ). which is insufficient to observe any of the effects 
discussed in this paper.) These planets have precession 
periods of around 100 years so that the argument of peri- 
apse of these planets changes by —10° during the course 
of Kepler observations. Though none of these planets lie 
in the Kepler field, they are all good candidates for ob- 
serving apsidal precession though precision photometry. 
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WASP-12b WASP-4b CoRoT-lb OGLE-TR-56b TrES-3b HAT-P-7b TrES-2b WASP-14b 



Fig. 6. — Eccentricities Needed to Detect Interior Properties from Apsidal Precession. The best-known planets for detecting 
k 2p precession are analogs to the hot Jupiters WASP-12b, WASP-4b, CoRoT-lb, OGLE-TR-56b, TrES-3b, HAT-P-7b, TrES-2b, and 
WASP- 14b. Assuming that analogs to these planets exist in the Kepler field around a V=14 magnitude star, the above graph shows the 
eccentricities required to detect k2 p . Black symbols correspond to calculations with u) = 0° and gray symbols correspond tou = 90°; in 
both cases, 6 = 0. Apsidal precession is much easier to detect for larger eccentricities so increasing e decreases the detectable k2 p . Using our 
transit light curve model, we found that threshold ki v values followed a power law k2 p oc e" 1 (for low eccentricities), which is consistent 
with the analytical estimates that -S- <x eu> <x ek2 p (see Section 13.21 1. Interpolating (and sometimes extrapolating) on this power law 
relationship, the graph identified the eccentricities required of these analog planets to detect precession due to a "typical" planetary interior 
of &2p = 0.3 (triangles). For example, when e = 0.00026 and lu = 0°, the apsidal precession due to an analog of WASP-12b should be just 
detectable by Kepler. A higher eccentricity (shown in Table 1) would be needed to measure &2p with sufficient accuracy (0.1) to distinguish 
between a massive core and a core-less model (circles). Systematic errors are expected to become important once the measurement error 
on &2p reaches as low as 0.01 (squares). If any of the very hot Jupiters discovered by Kepler have comparable eccentricities, the long-term 
high-precision photometry would allow for a powerful probe into their interior structure. HAT-P-7b and TrES-2b are known to lie in the 
Kepler observing field, but the values above are not corrected for improved photometric accuracy obtainable on these bright stars. Note 
that the eccentricities shown above and in Table 1 are computed for jL = 1; 3-cr measurements require eccentricities 3 times as high. 



WASP- 14b is more massive and has a larger semi- 
major axis (0.035 instead of 0.025) which is enough to 
significantly reduce the detectability of apsidal preces- 
sion which only proceeds at 0.1° per year. Unlike the 
previously mentioned planets, WASP - 14b has a know n 
non-zero eccentricity of 0.091 ± 0.003 (jJoshi et al.ll2009T i. 
Thus, the amplitude of transit timing variations is known 
to be very large (~97 minutes), but with a ~3400 year 
precession period. 

CoRoT-7b is a very hot super-Earth and has the short- 
est known or bital period (exce pting the ultra-short pe- 
riod planets of lSahu et al.ll2006l ). We included this planet 
in our analysis to get a feel for the plausibility of detect- 
ing the interior structure of terrestrial extra-solar plan- 
ets. The small radius reduces the planetary contribu- 
tion to apsidal precession (Figure 1) and significantly 
reduces the photometric signal. We note here that in 
bodies where material strength (rigidity) is more im- 
portant than self-gravity, fc2 P is no longer directly re- 
lated to internal density distribution. The correction 
factor is typically small fo r bodies larger than the Earth 
(|Murrav fc Dermottl[l999l) . 

XO-3b is a super-massive eccentric planet that is not 
in the Kepler field. Even so, it is interesting to note 
that, using the kno wn eccentricity e = 0.2884 ± 0.0035 
l|Winn et al.ll2009bT ) and accounting for the brightness of 
the host star (V=9.8), the Kepler thr eshold ki v is re- 
duced to o nly 0.54. As pointed out by Uordan h Bakosl 
(|2008f) and iPal fc Kocsisl (J2008D . XO-3b is a good can- 
didate for observing apsidal precession within the next 



decade or so. Furthermore, as discus sed below, the non - 
zero obliquity of the stellar spin axis (|Winn et al.ll2009bt ) 
may also result in an observable signal due to nodal pre- 
cession. 

HAT-P-7b and HAT-P-llb are orbiting bright stars 
in the Kepler field. The latter is an eccentric hot Nep- 
tune with a relatively large semi-major axis resulting in 
no eminently detectable apsidal precession. HAT-P-7b, 
on the other hand, is a good candidate for detecting 
apsidal precession. It is probably one of the brightest 
hot Jupiters in the Kepler field, orbiting a V=10.5 star. 
The system brightness improves the expected photomet- 
ric accuracy from 1000 ppm/min to 200 ppm/min, im- 
plying that an eccentricity of only 0.014 is needed to 
detect apsidal precession (threshold /c2p=0.3). IPal et all 
(|2009h report a best-fit eccentricity of 0.003 ± 0.012, in- 
dicating that the necessary eccentricity cannot be ruled 
out. Furthermore, this planet has transiting data extend- 
ing back to 2004 and was observed by NA SA's EPOXI 
Mission in 2008 (|Christiansen et all 120091: D. Deming, 
pers. comm.). This additional baseline, though sparsely 
sampled, may provide the additional leverage needed to 
detect apsidal precession if the eccentricity is non-zero. 
Note, however, that detecting changes in transit shapes 
is more difficult when the observations are made with a 
variety of telescopes because transit shapes depend on 
the observing filter used, due to wavelength-dependent 
limb darkening. 

TrES-2b is similar to HAT-P-7b in that it also lies in 
the Kepler field, has observations dating to 2005, and 
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TABLE 1 

Extra-Solar System Parameters and Results 



Planet Analog 


M. 


R* 


M p 


Rp 


a 


i b 


(JJtot 


e (Threshold fc 2p =0.1) c 


Threshold P d 


Threshold Q, d 


Ref 




M Q 


«© 


Mj 


Rj a 


AU 


ppm 


°/yr 


w = 0° 


w = 90 o 




ms/yr 










WASP-12b 


1.35 


1.57 


1.41 


1.79 


0.0229 


1830 


19.9 


0.0008 


0.0004 




0.95 






92700 


1 


CoRoT-lb 


0.95 


1.11 


1.03 


1.55 


0.0245 


314 


4.96 


0.0028 


0.0085 




0.93 






12500 


2,3 


WASP-4b 


0.92 


0.91 


1.24 


1.36 


0.0234 


109 


2.91 


0.0047 


0.0394 




0.68 






9900 


4 


TrES-3b 


0.93 


0.83 


1.91 


1.34 


0.0228 


106 


2.04 


0.0062 


0.0614 




0.53 






13700 


5 


OGLE-TR-56b 


1.17 


1.32 


1.29 


1.30 


0.0236 


451 


3.00 


0.0077 


0.0096 




1.36 






24700 


6 


HAT-P-7 b 


1.47 


1.84 


1.77 


1.36 


0.0377 


176 


0.25 


0.2085 


0.3146 




6.73 






2800 


7 


TrES-2 b 


0.98 


1.00 


1.19 


1.22 


0.0367 


18 


0.13 


0.2102 


e 




2.94 






350 


8 


WASP- 14b 


1.21 


1.31 


7.34 


1.28 


0.0360 


144 


0.09 


0.8352° 


c 




3.92 






5400 


9 


XO-3 b 


1.21 


1.37 


11.8 


1.22 


0.0454 


46 


0.04 


c 


c 




8.00 






1700 


10 


HAT-P-llb 


0.81 


0.75 


0.081 


0.42 


0.0530 


0.2 


0.01 


c 


c 




29.2 






0.1 


11 


CoRoT-7b 


0.91 


1.02 


0.028 


0.16 


0.0170 


8 


0.29 


c 


c 




16.8 






80 
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system parameters were used to estimate the 


detectability of aps 
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for 


these very 


hot Jupiter systems 
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For all systems 
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limb darkening paramctci 
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0.4 (appropriate 


for Kepler's bandpass) were used 


IMandel & 


Ago! 2002). 


For reference, the 



measured eccentricity of WASP-12b, WASP-14b, HAT-P-llb, and XO-3b arc 0.049 ± 0.015, 0.091 ± 0.003, 0.198 ± 0.046, and 0.2884 ± 0.0035 
respectively. Other planets have unmeasured eccentricities or eccentricity upper limits of <0.05. A discussion of these results is provided in 
Section [3~4| 

a We use Rj = 71492 km, the equatorial radius at 1 bar. The estimated depth of the secondary transit in Kepler's bandpass (sec Section 
|3.1|| . c The eccentricity required (at two different values of cj) so that a Ai2 p difference of 0.1 has an effective signal-to-noise of 1 in all of Kepler 
data for a V— 14 star, corresponding to a photometric accuracy of 1000 ppm/min. If analogs to these planets were found by Kepler with the given 
eccentricities, the internal density distribution would be measured well enough to detect the presence of a large core (see Section |3.2|| . These 
values correspond to the circles in Figurc[6] These results are for central transits (for b > 0, sec Figurc[7}- The value of the change in period, 
P, that can be detected with a signal-to-noisc of 1 in all of Kepler data for a V— 14 star (see Section |4.2|| . The value of threshold Q* is an 
estimate of the maximum value of the stellar tidal dissipation parameter, Q*, assuming that the period decay is due entirely to tidal evolution 
of the planet. Lower values of Q* are detectable by Kepler. Stars are thought to have time-averaged Q* values around 10000, though this value 
is highly uncertain and could be much higher for individual stars. 6 Even with the precision of Kepler, apsidal precession for these planets is 
undetectable. The extrapolation used to compute eccentricities at specific values of threshold k^ p assumes the inverse relationship discussed in 
the text &2p oc e~ , which is only true for low eccentricities. This ultra-short period low-mass planet was recently announced by the CoRoT 
team, but has not been published in a pccr-rcvicwcd journal. We take the parameters from J. Schneider's Extra-solar Planets Encyclopedia and 
use the mass-radius relation for terrestrial super-Earths of ISotin et al.| lj2 00 7f) to estimate the mass as ~9 Earth masses (rather than using the 
quoted upper limit of 17 Earth masses). 



was observed by NASA's EPOXI Mission. TrES-2b is 
somewhat fainter than HAT-P-7b (V=11.4), and, cor- 
recting for the system brightness, an eccentricity of 0.021 
would result in detectable apsidal precession (threshold 
fc2p=0.3). Observations of the secondary eclipse show 
no detectable deviation s of the orbit from circularity 
(|Q'Donovan et all [20091. Even so, the light curve of 
this planet is quite sensitive to perturbations as it has 
a quite high impact parameter b = 0.854. Accounting 
for this impact parameter does not significantly change 
the required eccentricity. 

We conclude that Kepler may detect the cores of very 
hot Jupiters and probe their interior structure though 
their evolving transit light curve if eccentricities are 
above ^0.003. As future observations provide longer 
baselines for these observations, the sensitivity to inte- 
rior structure measurements will increase dramatically, 
significantly lowering the eccentricity needed to observe 
apsidal precession. 

In cases where apsidal precession is not observed, the 
data can set strong upper limits on planetary eccentric- 
ities. An upper limit on the eccentricity can be inferred 
by assuming that the planet has the minimal physically- 
plausible value of k 2p ~ 0.1. Null detections of apsidal 
motion should therefore provide upper limits on eccen- 
tricity comparable to the values shown in Table 1 (also 
shown by circles in Figure [S]). Such strong eccentricity 
constraints are valuable for improving our understanding 
of these close-in planets. 



4. POTENTIAL CONFUSION OF THE APSIDAL 
PRECESSION SIGNAL 

In the above, we have assumed that measuring Co is 
tantamount to measuring k 2p . This is justified by not- 
ing that the conversion Co to k 2p involves only factors 
that are very well characterized. In Section 2 and Figure 
1, we showed that k 2p is usually the dominant source of 
apsidal precession. The effects of k 2 * and general relativ- 
ity are well-understood and can typically be subtracted 
away without introducing serious uncertainty, even when 
they dominate the apsidal precession rate. From Equa- 
tion [51 converting the remaining Co p to k 2p requires only 

knowing j^-, e, — -, and n. The latter two are very ac- 
curat ely measured with even a few transi t light curves 
(e.g., iTorreset al.l 120081 : lSouthworthll2008f ). The eccen- 
tricity only enters the equation through the /2(e) and 
32(e) eccentricity functions (Equations [71 and II ip . and 
Kepler observations of secondary eclipse are sufficiently 
accurate to remove any systematic error due to these 
terms unless the eccentricity is large (e > 0.3). Deter- 
mining the mass ratio requires well-sampled radial ve- 
locity observations. The systems detected by Kepler are 
bright enough to get good mass measurements, especially 
since very hot Jupiters have large radial velocity ampli- 
tudes (K ~ 200 m/s). 8 The anticipated error in the 

8 Other than determining the mass ratio and constraining the ec- 
centricity, radial velocity information is thought to have a negligible 
contribution in constraining apsidal precession unless a serious ob- 
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Fig. 7. — Effect of Impact Parameter on Precession Sig- 
nal. The detectability of apsidal precession depends on the im- 
pact parameter (b) of the orbital track across the star. For u = 0° 
(solid), the signal of primary transits are most important, with 
transit shaping playing the largest role. (See Figure [3]) However, 
the strength of transit shaping is a function of i mpact parameter 
with the minimum effect analytically estimated by Jordan & Bakos 
J2008h and lPal fc Kocsisl <20(jcl) to be b = X/y/2 (vertical solid line). 
Using a full photometric model, we see the expected decrease in 
the shaping signal (i.e. requiring a larger k,2 p to reach j, = 1). 
Note that the signal is nearly maximal, with small threshold k2 p 
values, for a large range of impact parameters. When w = 90° 
(dotted), the effect of primary transits are minimal and the offset 
in secondary transits become the determining factor. (See Figure 
U) At high impact parameters secondary eclipses are grazing, re- 
ducing the observable signal. We also show the threshold k2 p for 
an orientation of w = 45°, which lies, as expected, between the 
two extremes. The values of threshold k2 p shown are for an V=14 
CoRoT-lb analog in the Kepler field with an eccentricity of 0.003. 



mass ratio is a few percent ((Torres et al.l l2008h . In all, 
we estimate that, converting from Co to fc 2p leads to a 
typical systematic error on fc 2p of around ~.01. This is 
a relatively small systematic effect in comparison to the 
potential range (~0.5) of fc 2p values. For reference, the 
eccentricity required to reach a threshold fc 2p of 0.01 is 
shown in Figure [6] by squares. 

Another way to introduce systematic errors on the 
measurement of /e 2p is to misinterpret similar transit light 
curve variations. To ensure that the method outlined in 
this paper truly probes the interiors of extra-solar plan- 
ets, we consider in this section whether the transit light 
curve resulting from apsidal precession can be confused 
with any other common circumstances. Although a very 
specific combination of parameters is required for any 
particular phenomenon to successfully mimic a signal due 
to fc 2p , the below effects should be reconsidered when ac- 
tual data is available. 

4.1. Testing the Effect of Obliquity 

If either the star or planet has a non-zero obliquity, 
the orbital plane will no longer be fixed as a result of 
nodal precession. The obliquities of very hot Jupiters 

servational campaign can measure the radial velo city period (indc- 
pendently of transits) to sub -second accuracies. (Hcvl & Gladmanl 
2007: Jordan & Bakos 2008). 



rapidly (< 1 MYr) decay to a Cassini state, and re- 
cent work has s hown that these planet s are likely in 
Cassini state 1 flW inn & Holman 2005; iLevrard et al.l 
120071: iFabrvckv et al.l l2007t). Using a model b a sed o n 
the equations of lEggleton fe Kiseleva-Eggletonl (2001), 
we found that Cassini obliquities of very hot Jupiters are 
indeed negligible (a p < 0.01°). Though tidal damping 
of the stellar obliquity occurs on far longer timescales, 
several measurements of the projected stellar obliquity 
through the Rossiter-McLaughlin effect indicates that 
planet-hostin g stars generally have low obliquities < 10° 
like the Sun (|Fabrvckv fc Winnll2009D . Hence, the gen- 
eral expectation is that both the star and planet will have 
rather low, but potentially non-zero obliquities. 

Understanding the specific orbital evolution resulting 
from non-zero obliquities is more complicated than the 
simple prescription for apsidal precession. To correctly 
account for non-Kepleri an effects, we wr ote a direct inte- 
grator, following iMardling fc Linl (|2002f ). that calculates 
the Cartesian trajectory (and the direction of the spin 
axes) of a star-planet system including general relativ- 
ity and the effects of quadrupolar distortion. This inte- 
gr ator reproduces the orbit-averaged analytic equations 
of Mardling & Lin ( 200 2|), which are t h e same as thos e 
in lEggleton fc Kiseleva-Eggletonl (|2001h . lSterne] (|1939al) . 
and elsewhere. 9 We did not include the effects of tidal 
forces or additional planets which are not relevant to our 
problem. 

Using this direct integrator, we investigated the effect 
of non-zero obliquities on the transit times, durations, 
and impact parameters. Integration of several cases with 
varying stellar and planetary obliquities showed that the 
largest effect on the photometry was due to changes in 
the impact paramete r, as expected for an orbit with 
changing orientation (|Miralda-Escudel r2002f) . However, 
even for large stellar obliquities (~ 45°) the transit light 
curve variations due to obliquity are generally small rela- 
tive to the effects of purely apsidal precession, even with 
low eccentricities. One reason for this is that the tidal 
bulge, which does not contribute to nodal precession, is 
> 15 times more important than the rotational bulge. 
As with apsidal precession, the planetary contribution to 
orbital variations is much stronger than the stellar con- 
tribution (for equal obliquities). Unless the planetary 
obliquity is unexpectedly large (> 0.5°), the obliquity- 
induced nodal precession should have only a minor effect 
on the transit light curve. 



4.2. Transit Timing due to Orbital Decay 

Orbi tal decay genera tes a small secular trend in transit 
times. ISasselovl (|2003f) proposed the detectability of the 
expected ~1 ms/yr period change due to semi-major axis 
decay of OGLE-TR-56b. The transit timing anomaly 
due solely to orbital decay (or growth) is the result of 
constantly accumulating changes in the period: 



T N 



T + NP ohs + -N'dP 



(23) 



9 This involved mi nor modification s to th e "direct integrator" 
equations 3 and 5 in [Mardling & Lin (2002). In Equation 3, the 
coefficient 12 should be a 6 (R. Mardling, pers. comm.) and Equa- 
tion 5 was replaced with the nearly equivalent equation from Soffcl 
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where 5P = PP is the change in the period during 
one orbit and N is the number of transits after the 
initial transit. Equation [521 can be derived by not- 
ing that the transit times are basically the integral of 
the instantaneous period. As before, the transit tim- 
ing anomaly is composed of the quadratic deviation of 
Tjv from a straight line. The change i n period can be 
due to magnetic stellar breaking (e.g., iLee et al.l [20091 : 
iBarker fc Qgilvie 2009a), the Yarkovsky effect applied to 
planets (|Fabryckvil2008a( ) . and/or other effects. 

For planets orbiting an asynchronously rotating star, 
a major source of orbital decay is tidal evolution, which 
results in a s low change in semi-major axis, according to 
the formula ([Murray fc Dermottl[l999l) : 



. , , 3fo* M p 



(24) 



where sign(x) returns the sign of a: or if a; = and where 
Q* is the tid al quality parameter of the star, typically 
around 10 4 (|Dobbs-Dixon et al.ll200l . Though 5P due 
to tidal dissipation is only of order 3 micro-seconds, N 
grows by ~ 300 each year, reaching ^1000 during the 
duration of Kepler for very hot Jupiters. This implies a 
transit timing signal of about a few seconds. 

Calculating the total "signal-to-noise" of tidal evolu- 
tion, as was done for fo p , we find that reasonable values 
of Q* can be measured even for faint stars (V = 14; 1000 
ppm/min noise). For a circular orbit with the parameters 

of OGLE-TR-56b, the effective -£ reaches 1 when P is 
1.36 ms/yr (see Table 1), corresponding to Q» w 25000. 
This implies the detect ability of most o f the empirically- 
motivated estimates of Sas selovl (|2003| ) for the tidal de- 
cay of OGLE-TR-56b, which are estimated to be within 
an order of magn i tude o f 1 ms/yr. On the other hand, 
I Barker fc Ogilvia (|2009bf) estimate that the tidal damp- 
ing in F-stars like OGLE-TR-56 and WASP-12 may be 
very low, which may explain the survival of these short 
period planets. 

The estimates of the threshold values of P, shown in 
Table 1, include removing degeneracies in other param- 
eters, except apsidal precession of eccentric orbits, and 
assume that everything but P is known. Note that the 
transit light curve signal due to orbital decay is due en- 
tirely to transit timing; the change in a is far too small to 
observe in transit shaping. As the signal due to apsidal 
precession includes significant changes to the shapes of 
the transits, the signal due to fo p is qualitatively different 
than that of Q„. The shifting of secondary transits from 
precession also help in this regard, as outlined above. 
However, the primary transit timing signals can be simi- 
lar: quadratic transit timing anomalies with amplitudes 
of ~1 second. 

Kepler analogs of very hot Jupiters WASP-12b, 
OGLE-TR-56b, CoRoT-lb, WASP-4b, and TrES-3b 
could have detectable transit timing anomalies due to 
tidal decay, implying a direct measurement of the cur- 
rent value of Q* for specific stars (Table 1). This is an 
exciting possibility, providing the first direct measure- 
ments (or constraints) of the currently unknown details 
of tidal dissipation in a variety of individual stars. 10 We 

10 The vanishingly small effect of eccentricity decay is ~ -X- 



also note that interesting orbital decay of eclipsing binary 
systems seen by Kepler could also be detectable. 

4.3. Confusion Due to Other Planets 

Could the signal due to fo p be confused with addi- 
tional planets? In considering this issue, it should be 
noted that all known hot Jupiters (with a < 0.05 AU 
and M p > 0.5Mj up ) have no currently known additional 
companions. The apparent single nature of the se systems 
could very well be due to observational biases (|Fabrvckvi 
l2008bf) . However, even for stars that have been observed 
for many years with radial velocity (e.g. 51 Peg, HD 
209458), there appears to be a strong tendency towards 
hot Jupiters as the only close-in massive planets. 

Previous studies of transit timing variations fo- 
cus on the effects of a dditional plan etary per- 
turbers (e.g., iHolman fc Mu rray 2005; Agol et all 120051 : 
iFord fc Holmanl l2007t iNesvornv fc Morbidellil 120081) . 
These authors find that nearby massive planets or even 
low-mass planets in mean-motion resonances would cause 
strong transit timing variations that are easily dis- 
tinguishable from the comparatively long-period tim- 
ing anomalies due to fo p . Relatively distant com- 
panions or non-resonant low-mass planets, however, 
can ind uce a linear ap sidal precession signa l just 
like fo„ (iMiralda-Escudel 120021: iHevl fc Gladmanl 120071 : 
I Jordan fc Bakosll2008h . The precession rate induced by 
a perturbing body is a function of its mass and semi- 
major axis. The interior structure of very hot Jupiters 
causes apsidal precession as fast as a few degrees per 
year. To match this precession rate would require, for 
example, another Jupiter-mass planet at < 0.1AU or 
a solar-mass star at ~1 AU. Even perturbers an order 
of magnitude smaller than these would be readily de- 
tectable using radial velocity observations and/or high- 
frequency transit time variations. When restricted to 
planets that are undetectable by other means, adding 
the precession due to the unknown perturbing planet 
would lead to an insignificant overestimate of fo p for 
very hot Jupiters. 11 When observing transiting planets 
with larger semi-major axes (a > 0.05 AU), the strength 
of planetary induced apsidal precession is reduced to a 
level comparable t o apsidal precession fr om a low-mass 
perturbing planet ([Jordan fc Bak os 2008) and confusion 
may be possible in these cases. 

Since the transit timing signal for apsidal precession is 
similar to a sinusoid, another potential source of confu- 
sion would be light-trave l time o ffsets due to a distant 
orbiting companion (e.g. JDeeg e t al. 2008). The transit 
timing signal due to stellar motion about the barycenter 
can be distinguished from fo p precession 12 by consider- 
ing the changes in transit shapes and primary-secondary 
transit time offsets, which are not affected by distant 
companions. 

We conclude that transit timing effects from other 

smaller than apsidal precession, so that direct measurements of Q p 
from eccentricity decay are not feasible. 

11 Conversely, as a consequence of the fast precession of very hot 
Jupiters due to their (unknown) interiors, it will be very difficult 
to detect the presence of additional perturbing planets in these 
systems from apsidal precession alone. 

12 Transit time anomalies due to Q„ (Section l4.2l l. however, can 
be confused with barycenter light-travel time shifts due to a distant 
planet that may be undetectable in radial velocities. 
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planets can be readily distinguished from the effects of 
apsidal precession. To address the issue of the tran- 
sit shaping signal due to additional planets, we wrote 
a simple three-body integrator (similar to the integra- 
tor mentioned above) to investigate the kinds of transit 
light curve signals created by additional planets. For the 
vast majority of additional planet parameters, the tran- 
sit timing deviations always carry far more signal than 
the minor deviations due to changes in the angular ve- 
locity 13 (ftr) or impact parameter (b), which together 
determine the transit shape as described in Section 13.31 
above. Generally, it is much easier to delay a transit by 
5 seconds than it is to shift the apparent transit plane 
by an appreciable amount. 

However, when the perturbing planet is on a plane 
highly- inclined to the transiting planet, changes in the 
transit shape can become detectable, even while the tran- 
sit timing variations are negligible. For example, a per- 
turbing planet of mass 10 -5 M* at 0.1 AU with a mutual 
inclination of 45° caused very hot Jupiter transit dura- 
tions to change by ~1 second/year. This kind of signal 
is the result of nodal precession induced by the perturb- 
ing pl anet, as originally pointed out by iMiralda-Escudel 
(2002). In our investigation, we found that the three- 
body nodal precession alters the impact parameter (b) 
but does not significantly affect the orbital angular ve- 
locity (ftr)- Conversely, the transit shaping signal due to 
&2p is generally produced by changes in both b and ft r , 
but at near-central transits, the effect of changing orbital 
velocity is dominant (see Section 13. 3p . In high-precision 
transit light curves, both the angular velocity and the 
impact parameter can be independently measured and 
hence the signals of apsidal and nodal precession are usu- 
ally distinct for all but the most grazing transits. 

Given the uniqueness of the apsidal precession signal 
induced by the planet's interior, it appears that if ad- 
ditional planets are not detectable in radial velocities, 
transit timing variations, or nodal precession, then they 
will not contribute to a misinterpretation of an inferred 
value of /c2 P for very hot Jupiters. Nevertheless, future 
measurements of ki v should check that these issues are 
unimportant within the context of the specific system 
being studied. 

Finally, we estimate that moons or rings with enough 
mass to bias an inferred kip would cause other read- 
ily detectable phot ometric a nomalie s (e.g., plan et-moon 
barycentric motion ISartoretti fc Schneiderlll999T ). In ad- 
dition, extra-solar moons with any significant mass are 
tidally unstable, especi ally around very hot Jupiters 
(|Barnes fe O'Brienll2002T) . 

5. OTHER METHODS FOR DETERMINING K 2 p 

5.1. Secular Evolution of a Two Planet System 

Measuring k? for a n ext ra-solar planet was suggested 
by IWu fc Goldreichl (|2002f) for the inner planet of HD 
83443. Unfortunately, later analyses have indicated 
that the supposed second planet in this system was 
actually an artefac t of the sparse radial velocity data 
(|Mavor et al.1 [20041. Nevertheless, this technique could 
be applied to other eccentric planetary systems with sim- 

13 The angular velocity is directly related to the star-planet sep- 
aration through conservation of angular momentum: r/ 2 . 



ilar conditions (|Mardling|l200l . IWu fe Goldreichl ([2001 ) 
showed that in a regime of significant tidal circularization 
and excitation from an additional planet, the ratio of ec- 
centricities depends on the precessio n rate which is domi- 
nated by &2p as shown above (see also lAdams fc Laughlinl 
I2006L who do not include precession due to the planetary 
quadrupole) . In theory, the current orbital state of such 
multi-planet systems gives an indirect measurement of 
the apsidal precession rate. 

5.2. Direct Detection of Planetary Asphericity 

Another method for determining interior properties 
of transiting planets would be to directly measure the 
asphericity due to the rotational or tidal bulge in pri- 
mary transit photometry. The height of the rotational 
and tidal bulges are q r h,2R P and (ft/^-Rp, respectively, 
where q r and qt are the dimensionless small parameters 
defined in Equation [2] and hi is anothe r Love number 
which, for fluid bodies, is simply k^ + 1 (Sterne 1939a). 
These bulges cause the disk of the planet to be slightly 
elliptical, subtly modifying the p hotometric signal, as 
discu ssed for rotational bulg es by ISeager fc Huil (2002) 
andlBarnes fc Fortn cv (2003). However, as discussed by 
iBarnes fc Fortnevl (|2003h . in real systems with actual ob- 
servations, the size of the rotational bulge is very difficult 
to determine as it is highly correlated with stellar and or- 
bital parameters that are not known a priori, e.g. limb 
darkening coefficients. 

The tidal bulge, whose height is also set by hz p , does 
not suffer from some of the difficulties involved with 
measuring the rotational bulge. It has a known ori- 
entation (pointing towards the star) so there is no de- 
gener acy from an unknown obliquity (Ba rnes fc Fortnevl 
2003). (Note, however, that for hot Jupiters, the obliq- 
uities must be tidally evolved to nearly zero, so this isn't 
really a problem with the rotational bulge.) In addi- 
tion, the signal due to oblateness is only significant near 
ingress/egress, but the tidal bulge is continuously chang- 
ing orientation throughout the entire transit. Though 
the tidal bulge is typically three times larger than the 
rotational bulge (Equation 2), the projection of the tidal 
bulge that is seen during a transit is small, proportional 
to sm.9 where 9 is the angle between the planet posi- 
tion and the Earth's line of sight. For very hot Jupiters 
that have semi-major axes of only <6 stellar radii, sm.9 
during transit ingress/egress reaches > i so that the pro- 
jected tidal bulge is about half as large as the rotational 
bulge. The extra dimming due to the tidal bulges (and 
rotational bulges) is as high as 2 x 10~ 4 for some plan- 
ets that are expected to have tides over 2000 km high 
(e.g. WASP-12b, WASP-4b, Corot-lb, OGLE-TR-56b); 
this compares very favorably with the photometric accu- 
racy of binned Kepler data at about 10 ppm per minute. 
However, we expect that, as with the rotational bulge 
alone, the combined signal from the rotational and tidal 
bulge will be highly degenerate with the unknown limb- 
darkening coefficients, as the size of the projection of the 
tidal bulge also varies as the distance to the center of the 
star. 

We note that using multi-color photometry should 
significantly improve the prospects of detecting non- 
spherical planetary transits since it breaks most of t hese 
degeneracies. For example, iKnutson et al.l (|2007bl ) use 
HST to observe transits of HD 209458b in 10 wavelength 
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bands and measure the planetary radius with a rela- 
tive accuracy (between bands) of 0.003i?j, of the same 
level as the c hange in shape due to oblateness and the 
tidal bulge. iPont et all (|2007af) made a similar mea- 
surement for HD 189733b and reached even higher rela- 
tive accuracy. Combining such measurements with other 
data (e.g. primary transits in the infrared, where limb- 
darkening is much smaller) and a stellar photosphere 
mod el (to correctly correl ate limb darkening parameters 
as in lAgol fc Steffen 2007) could yield detections of plan- 
etary asphericity, especially in very hot Jupiters which 
have the largest bulges. 

One possible source of confusion in interpreting plane- 
tary asphericity is the thermally-induced pressure effects 
of an unevenly radiated surfac e. In non-synchronous 
planets, the thermal tidal bulge (|Arras fc Socratesll2009f) 
can shift the level of the photosphere by approximately 
an atmospheric scale height, about 10~ 2 or 10 -3 plane- 
tary radii (P. Arras, pers. comm.). The orientation of 
the thermal bulge is significantly different from the tidal 
or rotational bulges and should be distinguishable. Fur- 
thermore, very hot Jupiters should orbit synchronously, 
reducing the importance of this effect. Nevertheless, the 
effect of atmospheric phenomena on measurements of 
planetary asphericity should be considered. 

Though difficult to disentangle from other small pho- 
tometric effects, high-precision multi-color photometry 
may be another viable method for measuring ki p . This 
technique is complimentary to detecting ki p from apsi- 
dal precession since it does not require that the planet 
is eccentric, nor does it require a long time baseline. On 
some planets, the two methods could be used together as 
mutual confirmation of the planetary interior structure. 

6. CONCLUSIONS 

The planetary mass and radius are the only bulk phys- 
ical characteristics measured for extra-solar planets to 
date. In this paper, we find that the planetary Love 
number (fc2p, equivalent to J2) can also have an ob- 
servationally detectable signal (quadrupole-induced ap- 
sidal precession) which can provide a new and unique 
probe into the interiors of very hot Jupiters. In par- 
ticular, ki p is influenced by the size of a solid core and 
other internal properties. Core sizes can be used to in- 
fer the formati on and evolution of individual extra-solar 
planets Ce.g.. iDodson- Robinson fc Bodenheimerl 120091 : 
iHelled fe Schubertll2009D ~~ 

The presence of a nearby massive star creates a large 
tidal potential on these planets, raising significant tidal 
bulges which then induce non-Keplerian effects on the 
star-planet orbit itself. The resulting apsidal precession 
accounts for ~95% of the total apsidal precession in the 
best cases (Figure [T]). Hence, we find that the internal 
density distribution, characterized by ki pi has a large and 
clear signal, not to be confused with any other parame- 
ters or phenomena. We urge those modeling the interior 
structures of extra-solar planets to tabulate the values of 
ki v for their various models. 

Encouraged by this result, we calculated full photo- 
metric light-curves like those expected from the Kepler 
mission to determine the realistic observability of the in- 
terior signal. We estimate that Kepler should be able to 
distinguish between interiors with and without massive 
cores (A^2 P — 0.1) for very hot Jupiters with eccentrici- 



ties around e ~ 0.003 (Figure[B]). Eccentricities this high 
may occur for some of the very hot Jupiters expected to 
be found by Kepler, though these planets usually have 
highly damped eccentricities. Much stronger constraints 
on apsidal precession can be obtained by combining Ke- 
pler photometry with precise secondary transits observed 
in the infrared. In cases where apsidal precession is not 
observed, the data can set strong upper limits on plane- 
tary eccentricities. 

In analyzing Kepler's photometric signal of apsidal 
precession, we find that transit timing variations are an 
almost negligible source of signal, though transit tim- 
ing has been the focus of many observational and the- 
oretical papers to date. The effect of "transit shaping" 
has ~30 times the photometric sign al of transit timing 
for apsidal precession (see Figure [H iPal fe Kocsisll2008t 
I Jordan fc Bakosl 1200 8)). At orientations where transit 
timing and shaping are weakest, the changing offset be- 
tween primary and secondary transit times can be used 
to measure fc 2p (Figure [4} . It may also be possible to 
measure ki p from high-precision multi-color photometry 
by directly detecting the planetary asphericity in transit. 
Such a measurement does not require a long baseline or 
an eccentric orbit. 

Very hot Jupiters are also excellent candidates for de- 
tecting tidal semi-major axis decay, where we find that 
relatively small period changes of P ~ 1 ms/yr should be 
detectable. This could constitute the first measurements 
(or constraints) on tidal Q* for a variety of individual 
stars. We note that Kepler measurements of transit tim- 
ing and shaping for eclipsing binaries should also provide 
powerful constraints on stellar interiors through apsidal 
motion and binary orbital decay (due to tides, if the com- 
ponents are asynchronous). 

Accurately measuring the interior structure of distant 
extra-solar planets seems too good to be true. Nev- 
ertheless, the exquisite precision, constant monitoring, 
and 3.5-year baseline of the Kepler mission combined 
with the high sensitivity of transit light curves to small 
changes in the star-planet orbit make this measurement 
plausible. 

Our focus on Kepler data should not be interpreted to 
mean that other observations will be incapable of mea- 
suring k2 P . In fact, the opposite is true since the size of 
the apsidal precession signal increases dramatically with 
a longer baseline. Combining Kepler measurements with 
future ground and space based observations can create 
a powerful tool for measuring ki p . In the far future, 
many planets will have measured apsidal precession rates 
(like eclipsing binary systems have now) and inferred ki p 
values. Incorporating these measurements into interior 
models holds promise for greater understanding of all 
extra-solar planets. 



We thank Dave Stevenson, Mike Brown, Greg Laugh- 
fin, Oded Aharonson, Thomas Beatty, Phil Arras, Rose- 
mary Mardling, Re'em Sari, Alejandro Soto, Ian McEwen 
and Chris Lee for help and useful discussions. We espe- 
cially thank the referee, Dan Fabrycky, for helpful sug- 
gestions and discussions. DR is grateful for the support 
of the Moore Foundation. ASW is grateful for support 
from the National Science Foundation. This research has 
made use of NASA's Astrophysics Data System. 



IS 



Ragozzine & Wolf 



REFERENCES 



Adams , F. C., & Laughlin, G. 2006, ApJ, 649, 1004, 

|arXiv:astro-ph /0606349 
Agol, E., Cowan, JN. B., Bushong, J., Knutson, H., Charbonncau, 

D., Deming, D., & Steffen, J. H. 2009, in IAU Symposium, Vol. 

253, IAU Symposium, 209-215 
Agol, E., Steffen, J., Sari, R., fc Clarkson, W. 2005, MNRAS, 359, 

567, |arXiv:astro-ph/0412032| 
A gol, E., fc Steffen, J. H. 2007, MNRAS, 374, 941, 

arXiv:astro-ph/0610159 
Alonso, R., Aigrain, S., Pont, F., Mazeh, T., & The CoRoT 

Exoplanet Science Team. 2009, in IAU Symposium, Vol. 253, 

IAU Symposium, 91-96 
Arras, P., & Socrates, A. 2009, ArXiv e-prints, 0901.0735 
Bakos, G. A. et al. 2009, ArXiv e-prints, 0901.0282 
Baraffe, I., Chabrier, G., & Barman, T. 2008, 0802.1810 
Barge, P. et al. 2008, 0803.3202 
Barker, A. J., & Ogilvie, G. I. 2009a, ArXiv e-prints, 0902.4554 

. 2009b, ArXiv e-prints, 0902.4563 

Barne s, J. W., fc Fort ney, J. J. 2003, ApJ, 588, 545, 

arXiv:astro-ph/0301156 
Barne s7~J~ W., fc O 'Brien, D. P. 2002, ApJ, 575, 1087, 

|arXiv:ast ro-ph/0205035 
Batten, A. H. 1973, Binary and multiple systems of stars (Binary 

and multiple systems of stars / by Alan H. Batten. Oxford 

; New York : Pergamon Press, [1973] (International series of 

monographs in natural philosophy ; v. 51)) 
Bean, J. L. 2009, ArXiv e-prints, 0903.1845 

Beatty, T. G., & Gaudi, B. S. 2008, ArXiv e-prints, 804, 0804.1150 
Bodenheimer, P., Lin, D. N. C, & Mardling, R. A. 2001, ApJ, 548, 

466 
Borucki, W. J. et al. 2003, in Presented at the Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference, Vol. 

4854, Future EUV/UV and Visible Space Astrophysics Missions 

and Instrumentation. Edited by J. Chris Blades, Oswald H. W. 

Siegmund. Proceedings of the SPIE, Volume 4854, pp. 129-140 

(2003)., ed. J. C. Blades & O. H. W. Siegmund, 129-140 
Burrows, A., Hubeny, I., Budaj, J., & Hubbard, W. B. 2007, ApJ, 

661, 502, arXiv:astro-ph/0612703 
Burrows, A., Ibgui, L., & Hubeny, I. 2008, ArXiv e-prints, 803, 

0803.2523 
Charbonneau, D. 2003, in Astronomical Society of the Pacific 

Conference Series, Vol. 294, Scientific Frontiers in Research on 

Extrasolar Planets, ed. D. Deming & S. Seager, 449—456 
Charbonne au, D. et al. 2005, ApJ, 626, 523, 

arXiv:astro-ph/0503457 
Charbonneau, P., Bro wn, T. M., Latham, D. W., & Mayor, M. 

2000, ApJ, 529, L45, [arXiv:astro-ph/9911436| 
Chatterjee, S., Ford, E. B., fc Rasio, Jj\ A. 2007, ArXiv 

Astrophysics e-prints, astro-ph/0703166 
Christian, D. J. et al. 2008, ArXiv e-prints, 806, 0806.1482 
Christiansen, J. L. et al. 2009, in IAU Symposium, Vol. 253, IAU 

Symposium, 301—307 
Claret, A. 1995, A&AS, 109, 441 
Cowling, T. G. 1938, MNRAS, 98, 734 
Deeg, H. J., Ocana, B., Kozhevnikov, V. P., Charbonneau, D., 

O'D onovan, F. T., & Doyle, L. R. 2008, A&A, 480, 563, 

arXiv:0801.2186 
Dobbs -Dixon, I., Lin, D. N. C. fc Mardling, R. A. 2004, ApJ, 610, 

464, [arXiv:astro-ph/0408191| 
Dodson- Robinson, S. E., & Bodenheimer, P. 2009, ArXiv e-prints, 

0901.0582 
Egglet on, P. P., fc Kisele va-Eggleton, L. 2001, ApJ, 562, 1012, 

arXiv:astro-ph/0104126 
Fabrycky, D. 2008a, ApJ, 677, L117, 0803.1839 
Fabry cky, P., fc Tremaine, S. 2007, ApJ, 669, 1298, 

arXiv:0705.4285 
Fabrycky, D. C. 2008b, ArXiv e-prints, 806, 0806.4314 
Fabryc ky, D. C, J ohnson, E. T., & Goodman, J. 2007, ApJ, 665, 

754, ^rXiv:astro-ph/0703418| 
Fabrycky, D. C, & Winn, J.N. 2009, ArXiv e-prints, 0902.0737 
Ferraz-Mello, S., Rodriguez, A., & Hussmann, H. 20 08, 0712.1156 
Ford, E. B., & Holman, M. J. 2007, ApJ, 664, L51, arXiv:0705.0356 
Ford, E. B., Lystad, V., & Rasio, F. A. 2005, Nature, 131, 873, 

arXiv:astro-ph/0502441 



Gillon, M., Anderson, D. R., Demory, B. ., Wilson, D. M., Hellier, 

C, Queloz, D., & Waelkens, C. 2008, ArXiv e-prints, 806, 

0806.4911 
Goldreich, P., & Soter, S. 1966, Icarus, 5, 375 
Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, 

Nature, 435, 466 
Goodman, J. 2009, ArXiv e-prints, 0901.3279 
Gu, P.-G., & Ogilvie, G. I. 2009, ArXiv c-prints, 0901.3401 
Guillot, T. 2005, Annual Review of Earth and Planetary Sciences, 

33, 493, arXiv:astro-ph/0502068 
Guillot, T., Santos, N. PontTF., Iro, N., Melo , C, & Ribas, I. 

2006, A&A, 453, L21, arXiv:astro-ph/0605751 
Hebb, L. et al. 2009, ApJ, 693, 1920, 0812.3240 
Helled, R., & Schubert, G. 2009, ArXiv e-prints, 0903.1997 
H eyl, J. S., fc Glad man, B. J. 2007, MNRAS, 377, 1511, 

arXiv:astro-ph/0610267 
Holman, M. J., fc Murray, N. W. 2005, Science, 307, 1288 
Holman, M. J. et al. 2007, ApJ, 664, 1185, 0704.2907 
Hood, B., Wood, K., Seager, S., & Collier Cameron, A. 2008, ArXiv 

e-prints, 807, 0807.1561 
Hut, P. 1981, A&A, 99, 126 

Iorio, L. 2006, New Astronomy, 11, 490, |arXiv:gr-qc/0 505107 
Jackson, B., Greenberg, R., & Barnes, R. 2008, 0801.0716 
Johns-Krull, C. M. et al. 2008, ApJ, 677, 657, 0712.4283 
Jordan, A., & Bakos, G. A. 2008, ArXiv e-prints, 806, 0806.0630 
Joshi, Y. C. et al. 2009, MNRAS, 392, 1532, 0806.1478 

— . 2008, ArXiv e-prints, 806, 0806.1478 
Kallrath, J., Milone, E. F., Kallrath, J., & Milone, E. F., cds. 1999, 

Eclipsing binary stars : modeling and analysis 
Kjeldsen, H, Bedding, T. R., & Christensen-Dalsgaard, J. 2008, 

ArXiv e-prints, 807, 0807.0508 
Knutson, H. A. et al. 2007a, Nature, 447, 183, arXiv:0705.0993 
Knutson, H. A., Charbonneau, D., Noye s, R. W, Brown, T. M., fc 

Gilliland, R. L. 2007b, ApJ, 655, 564, |arXiv:astro-p h/0603542 
Koch, D. et al. 2006, Ap&SS, 304, 391 
Kopal, Z. 1959, Close binary systems (The International 

Astrophysics Series, London: Chapman & Hall, 1959) 
Kopal, Z., ed. 1978, Astrophysics and Space Science Library, 

Vol. 68, Dynamics of Close Binary Systems 
Kozai, Y. 1959, AJ, 64, 367 

Lanza, A. F. et al. 2009, A&A, 493, 193, 0811.0461 
Laughlin, G., Marcy, G. W., Vogt, S. S., Fischer, D. A., & Butler, 

R. P. 2005, ApJ, 629, L121 
Lee, J. W., Kim, S.-L., Kim, C.-H, Koch, R. H., Lee, C.-U., Kim, 

H.-L, & Park, J.-H. 2009, AJ, 137, 3181, 0811.3807 
Levrard, B., Correia, A. C. M., Chabrier, C, Baraffe, I., Selsis, F., 

& Laskar, J. 2007, A&A, 462, L5, arXiv:astro-ph/0612044 
Levrard, B., Winisdoerffer, C, & Chabrier, G. 2009, ApJ, 692, L9, 

0901.2048 
Loeb, A. 2005, ApJ, 623, L45, arXiv:astro-ph/0501548 

— . 2008, ArXiv e-prints, 807, 0807.0835 
Lopez-Morales, M., & Seager, S. 2007, ApJ, 667, L191, 

arXiv:0708.0822 
Mand el, K., fc A gol, E. 2002, ApJ, 580, L171, 

arXiv:astro-ph/0210099 
Mardling, R. A. 2007, MNRAS, 382, 1768, arXiv:0706.0224 
Mardling, R. A., & Lin, D. N. C. 2002, ApJ, 573, 829 
Matsumura, S., Takeda, G., & Rasio, F. A. 2008, ApJ, 686, L29, 

0808.3724 
Mayor, M., Udry, S., Naef, D., Pepe, F Queloz, P., Santos, N. C, 

& Burnet, M. 2004, A&A, 415, 391, arXiv:astro-ph/0310316 
Miralda-Escude, J. 2002, ApJ, 564, 1019, arXiv:astro-ph/0104034 
Murray, C. D., & Dermott, S. F. 1999, Solar system dynamics 

(Solar system dynamics by Murray, C. D., 1999) 
Nesvorny, D., & Morbidelli, A. 2008, ApJ, 688, 636 
O'Donovan, F. T., Charbonneau, D., Harrington, J., Seager, S., 

Deming, D., & Knutson, H. A. 2009, in IAU Symposium, Vol. 

253, IAU Symposium, 536-539 
Pal, A., Bakos, G. A., Noyes, R. W., & Torres, G. 2009, in IAU 

Symposium, Vol. 253, IAU Symposium, 428-431 
Pal, A., & Kocsis, B. 2008, ArXiv e-prints, 806, 0806.0629 
Peale, S. J., Cassen, P., & Reynolds, R. T. 1979, Science, 203, 892 
Pont, F. et al. 2007a, A&A, 4 76, 1347, larXiv:0707. 19401 

— . 2007b, A&A, 465, 1069, |arXiv-astr o-ph/0610827 
Rafikov, R. R. 2008, ArXiv e-prints, 807, 0807.0008 
Rodriguez, A., & Ferraz-Mello, S. 2009, ArXiv e-prints, 0903.0763 



Probing Extra-solar Planet Interiors 



19 



Rowe, J. F. et al. 2007, ArXiv e-prints, 711, 0711.4111 
Russell, H. N. 1928, MNRAS, 88, 641 

Sahu, K. C. et al. 2006, Nature, 443, 534, |arXiv:astro -ph/0610098 
Sartoretti, P., & Schneider, J. 1999, AfeA S, 134, 553 
Sasselov, D. D. 2003, ApJ, 596, 1327, arXiv:astro-ph/0303403 
Scharf, C. A. 2007, ApJ, 661, 1218, arXiv:astro-ph/0702749 
Seager, S., & Hui, L. 2002, ApJ, 574, 1004, arXiv:astro-ph/0204225 
Sing, D. K., & Lopez-Morales, M. 2009, A&A, 493, L31, 0901.1876 
Soffel, M. H. 1989, Relativity in Astrometry, Celestial Mechanics 
and Geodesy (Relativity in Astrometry, Celestial Mechanics 
and Geodesy, XIV, 208 pp. 32 figs.. Springer- Verlag Berlin 
Heidelberg New York. Also Astronomy and Astrophysics 
Library) 
Sotin, C, Grasset, O., & Mocquet, A. 2007, Icarus, 191, 337 
Southworth, J. 2008, 0802.3764 
Sozzetti, A. et al. 2009, ApJ, 691, 1145, 0809.4589 
Sterne, T. E. 1939a, MNRAS, 99, 451 

— . 1939b, MNRAS, 99, 662 
Swain, M. R., Vasisht, G., & Tinetti, G. 2008, Nature, 452, 329 
Thommes, E. W., Bryden, G., Wu, Y., & Rasio, F. A. 2008, ApJ, 
675, 1538, arXiv:0706.1235 



Torres, G., Winn, J. N., & Holman, M. J. 2008, 0801.1841 
Winn, J. N., fc Holm an, M. J. 2005, ApJ, 628, L159, 

|arXiv:astro-ph /0506468 

Winn, J. N. et al. 2007, AJ, 134, 1707, arXiv:0707.1908 

Winn, J. N., Holman, M. J., Carter, J. A., Torres, G., Osip, D. J., 

& Beatty, T. 2009a, AJ, 137, 3826, 0901.4346 
Winn, J. N., Holman, M. J., Shporer, A., Fernandez, J., Mazeh, 

T., Latham, D. W., Charbonneau, D., & Everett, M. E. 2008, 

AJ, 136, 267, arXiv:0804.2479 
Winn, J. N. et al. 2009b, ArXiv e-prints, 0902.3461 

— . 2006, ApJ, 653, L69, arXiv:astro-ph/0609506 

— . 2005, ApJ, 631, 1215, arXiv:astro-ph/ 0504555 
Wu, Y., fc Goldreic h, V. 20TJ2; ApJ, 564, 1024, 

arXiv:astro-ph/0108499 
Zharkov, V. JN., &: 'Itubitsyn, V. P. 1978, Physics of planetary 

interiors (Astronomy and Astrophysics Series, Tucson: Pachart, 

1978) 



